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The study of metal-semiconductor junctions started in 
1874 with the work of Braun [l]. Although more than a century 
has elapsed since then, the understanding of the metal- 
semiconductor junction is still far from complete. Perhaps 
the main reason for this is the fact that the performance of 
semiconductor junctions is strongly dependent on processing, 
which in turn controls the condition of the metal semiconductor 
interface. For example, it is commonly found that depending 
upon the condition at the interface, the rectifying metal- 
semiconductor junctions (Schottky harrier diodes) made on 
low-doped silicon, GaAs etc. show either a temperature 
independent or a temperature dependent ideality factor. The 
study in the present thesis has been mainly directed towards 
the understanding and clarification of some of the ambiguities 
in those cases where the ideality factor is temperature 
dependent. Such ambiguities have been detected in the 
published literature concerning the assignment of the 
mechanism of current transport in some specific metal- 
semiconductor diodes. The clarification of these ambiguities 
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has been brought out in the present work by incorporating the 
fluctuation of geometrical and material parameters at the metal- 
semiconductor interface [2], 

A description of the ambiguities under consideration 
here would become easier if we first consider the commonly 
used method for assignment of the mechanism of current trans- 
port in ’homogeneous' Schottky barrier diodes; the term 
'homogeneous' that there is no fluctuation of . I 

geometrical and material parameters across the metal- 
semiconductor interface. The method of assignment of the 
mechanism of current transport for homogeneous Schottky- 
barrier diodes is well known and was developed by Saxena [ 3]. 
This method consists of obtaining the current voltage (l-V) 
characteristics in the forward region at different temperature 
(T), and thereby plotting nT versus T from the slope of the 
lnl-y plot, where n is the ideality factor of the diode. 
Different mechanisms of current transport are identified by 
their characteristic nT versus T plot given in Fig. 1 ,2.where 
the curves I, II and III represent the cases where the 
thermionic emission (TE) is the mechanism of ourrent transport 
with unity ideality factor, with ideality factor greater than 
unity and with T p effect respectively, CJurves I1T am# ■V 
correspond to thermionic field emission (TEE) and field 
emission (EE) respectively. 

This method of assignment of the mechanism of current 
transport has been used by many workers. Some [ 4 ] have found 
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that in the case of metal- (low doped) silicon diodes operating 
above room temperature the experimentally obtained nT is indepen- 
dent of temperature. Thus the assignment of transport is made 
to be by field emission. This however is contrary to the common 
experience [5,6] that Schottky barriers made on low-doped 
silicon and operating at temperatures higher than room temperature 
have thermionic emission as the mechanism of current transport. 

In the present work we have demonstrated that in the presence of 
fluctuations in geometrical and material parameters (transverse 
to the metal-semicondu ctor interface) there is a definite 
possibility of obtaining a temperature independent plot of nT 
even in some of those cases where thermionic emission is the 
dominant mechanism of current transport. Thus one concludes 
that the correspondence of field emission with temperature 
independent value of nT is not unambiguous. The repercussions 
of this conclusion are obvious - the nature of the nT versus 
T plot is not uniquely relatable to a specific mechanism of 
current transport. 

The role of fluctuations mentioned above in current 
transport has been incorporated in the present work in an 
intuitive fashion [ 2] . We have observed that the fluctuations 
of which is a characteristic energy characterising the 
interface distribution for a diode showing I effect (see 
Levine [7]), can play a dominant role in current transport. 
Furthermore, we have assumed that the different areas of 
Schottky barriers have different E o ’s and the distribution is 
Gaussian such that 
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where 0 Bi is the harrier height of the ith ptrtchf or an applied 
voltage V. Using equations (l) and (2) the I-Y characteristics 
and hence were evaluated at different temperatures.- 

Then using the relation, 


/ 3 In I" \-1 nkT 
( “*TT — > = “q~ 


(3) 


the nT versus T plot was obtained. Depending on the values of 
E o and a , it has been possible to get temperature independent 
nT for some temperature ranges, although the mechanism of 
current transport in equations (1) - (3) is only thermionic 
emission. 

Similar calculations have been carried out for Schottky 
barriers with a thin oxide layer between the metal and the 
semiconductor. Fluctuations in the interfacial oxide thickness 
have been incorporated in this calculation. It has been observed 
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that the fluctuations in oxide thickness do not appreciably 
change the nature of the nT versus T plot. This and the 
previous calculations considering the fluctuations in are 
the subject matter of the first and the second chapters of 
this thesis. 

An experiment has been designed to verify the 

theoretical model mentioned above. For this a large number 

of isolated small area (homogeneous) diodes have been 

prepared and the nT versus T plot of each has been obtained. 

These show a straight-line plot with an intercept T^, 

where there is a scatter in the values of T for different 

o 

diodes. This has been true even in the case of diodes made 
on the same wafer. We have selected three or four diodes ,»f 
different values of T q and connected them in parallel. Then 
the nT versus T plot of the combination has been obtained 
and it has been found that over certain range of temperature 
nT becomes independent of T. The condition under which these 
observations have been made and other details of the experiment 
are given in Chapter III. 

Having noted that the fluctuations in physical 
parameters could cause the ambiguities in the determination 
of the mechanism of current transport we have proceeded 
further to determine some phenomenological parameters which 
would correlate T q with geometrical and material parameters 
of the interface. One such parameter, S Q , which is the 
characteristic energy describing the assumed exponential 



distribution of interface states has already been introduced 
earlier. Regarding E , we note that the I q effect is not 
unique to an exponentially distributed surface states as 
predicted by Levine [7] but to. the' constancy of the slope of 
lnQ SCT^B temperatures and bias, as noted by 

Crowell [ S] . Here Q s ^ is the total depletion charge in the 
semiconductor and 0 B is the barrier height. Crowell C8] 
has also noted that the potential profile of metal-P (thin )-TST 
diodes for a certain doping profile in the P-layer is such 
that one can obtain a linear lnQ go -0 B plot leading to T q 
eff e rt. 

It is interesting to note here that a similar 
potential profile as in the case of M-P (thin)-N-diodes, is 
also obtained in the case of Heine’s model [9] of a Schottky 
barrier. It Is well known that Heine’s model starts on a very 
different premise from Bardeen’s model [ 10) of Schottky barri ec 
and primarily considers that the surface-states are due to the 
tunneling of metal electrons into the forbidden bandgap of the 
semiconductor* A direct consequence of Heine’s assumption is 
that in this case the surface state charges are not localised 
in space, This non-dooalization of space-charge makes the 
potential profile similar to that of M-P ( thin )-NT-diodes. 

Hence we attempted to correlate the I Q effect with the 
parameters of tunneling in Heine’s model as developed in 
Pellegrini’s caloulations[l 1 ], The details of this correla- 
tion are given in the. fourth chapter. 
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In the penultimate chapter, a novel method to determine 

the doping density of a semiconductor wafer using the dependence 

of 1 on the current at which T is measured has been described, 
o o 

This method requires the measurement of only the forward X-V 
characteristics of a Schottky barrier diodes. A concluding 
discussion of the various observations made in the present work 
is given in the last chapter. 
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CHAPTER 1 


INTROIUCTI OI 


The present study is aimed at the investigation of 
current transport in nonhomogeneous rectifying metal- 
semi conductor contacts. The nonhomogenities that have been 
considered are due to the fluctuation of geometrical and 
material parameters transverse to the metal-semiconductor 
interface. 

To put this work in its proper perspective a flow- 
diagram of major developments in the study of metal-semicon- 
ductor contacts has been given in Pig. 1.1. It is worth 
noting that the first reported work [1] in this area was 
done in 1874 by Braun, who made point-contacts to semicon- 
ductors like copper sulphide and iron sulphide. It is 
generally considered, for example see [2], that the next 
significant step in this study was the recognition by 
Schottky et al. [3] of a 'potential barrier' as the cause of 
the asymmetry in the current -volt age characteristics, 
developments from this point onwards were rather fast and 
various classifications regarding the order of development 
are possible. We have arbitrarily decided to classify 
these developments under two heads, (i) the 'study of 
current transport' and (ii) the 'study of barrier heights' j 
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it is obvious, however, that both of these studies are highly 
interrelated and a conceptual separation between them can 
only be artificial. To rectify this situation several cross- 
connecting links have been shown in Pig. 1.1. 

1 • 1 Current Transport 

The current transport in a metal-semiconductor contact 
can be either (i) due to tunneling of carriers across the 
potential barrier or (ii) due to the transport of carriers 
over the barrier. The tunneling mechanism was the first one 
to be suggested and then was rejected as ambiguous y it was 
later rehabilitated as a valid mechanism in some cases, as 
would become clear from the discussion below. 

1.1.1 Tr ans port by tunneling a cross t he barrier 

Wilson [4] in 1932, tried to explain the rectifying 
action in terms of quantum mechanical tunneling of electrons 
across the barrier. However, it was shown by Davydov [5] in 
1938 that for the Schottky barrier diodes made in that period, 
the direction of rectification was opposite to the direction 
of rectification predicted by the tunneling theory of Wilson. 
Thj.s observation of bavydov [5] led to tunneling being 
completely ruled out as a mode of current transport in 
Schottky barriers. The situation, however, changed when 
Padovani and Stratton [6] and later Crowell and Rideout *[7] 
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showed, that in the case of diodes made with highly doped 
semiconductors, the current transport is due to tunneling, 
thus giving rise to two additional modes of current transport, 
viz. Field Emission (FE) , wherein the current flow is 
dominated by carriers tunneling close to Fermi level and 
Thermionic Field Emission (TFE), where the current flow is 
due to temperature aided tunneling, away from the Fermi -level. 
In Fig. 1.1, the block representing the works of Padovani 
and Stratton [6] and Crowell and Rideout [7], is separated 
from that representing Wilson’s [5] work by a block repre- 
senting the work done on tunneling in Metal -Insulator Metal 
(M-I-M) structures. Though the study in MIM structure is 
not explicitly related to the study of Schottky barriers, the 
tunneling theory used to study the Schottky barrier with 
highly doped semiconductor was developed in the context of 
MIM devices only [8,9]. Having traced the development in . 
the current transport model for tunneling across the barrier, 
we move over to look into the development in the models 
describing the current transport by 'emission over the 
barrier ' . 

1.1.2 Transport by emission over the barrier 

As is shown in Fig. 1.1 there are two distinct models to 
explain current transport by emission over the barrier. The 
first model for transport of carriers over the barrier was 
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advanced by Schottky [10] and independently by Mott [11], 
when they pointed out that the current-voltage relationship 
can be explained assuming that the carriers passed over the 
barrier by the normal process of diffusion and drift. The 
next model of current transport was proposed by Bethe, in 
1942 [12], who assumed the transport of electrons over the 
barrier to be by the process of thermionic emission. The 
models proposed by Schottky and Mot t ' [ 10,11] , and Bethe [12] 
are respectively the ’Diffusion theory' and 'Thermionic 
emission (TE) theory'. A generalized transport model called 
thermionic emission diffusion (TED) theory was developed by 
Crowell and Sze in 1966 [13]. A good treatment of these 
development has been recently given by Rhoderick [2]. 

1.1.3 Methods to distinguish betw e en transport by tunneling 
and by "emission 

It is obvious that an experimental method which in a 
fabricated Schottky barrier diode can determine whether the 
current transport in. it is dominated by tunneling through the 
barrier or by emission over the barrier would be of great 
importance. A search for this technique in 1950' s and 1960's, 
culminated in the experiments suggested by Saxena In 1969 [14]. 
These experiments of Saxena are of central importance to our 
work and are discussed in detail in Sec. 1.3. Here we only 
want to reiterate that Saxena' s work provides a way to dis- 
tinguish between the TE, TPE and PE mechanisms of transport. 
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1 . 2 Barrier Height 

The first model to evaluate the harrier height was 
proposed by Mott in 1938 [11]. According to Mott the barrier 
<J> g seen from the metal would be given by the difference 
between the work-function $ of the met al and the electron 
affinity x s of the semiconductor. However, the experimental 
measurements [15] of brought to notice that there were 
metal-semiconductor systems in which Mott’s conclusions were 
not valid-, actually there were systems where the barrier 
height was almost independent of the work function of the 
metal used. This, in a way led to the work of Bardeen, Heine 
and Inkson as shown in Fig. 1.1, 

The first model to explain the weak dependence of <J> B 
ond recognised the importance of the existence of surface 
states on the semiconductor and was proposed by Bardeen in 
1947 [16]. The Bardeen model assumes the existence of a thin 
insulating layer between the metal and the semiconductor and 
identifies a neutral level 1 in the forbidden energy gap 
of the semiconductor, to characterise the surface states. 

It has been further demonstrated by Bardeen that in the case 
of large surface state densities, the semiconductor Fermi- 
level at the surface gets pinned to the neutral level M 0 > . 

Two of the important extensions of Bardeen's model are 
the study of Metal -Insulator -Semi conductor Structure (MIS) 

(by Cowley and Sze [17], and by Card and Rhoderick [18]) and 
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the exponential surface state model (by Levine [19]). The 
details of these two models are given in Sections 1.4 and 2.1 
respectively. The other interesting extension of Bardeen’s 
work, shown in Big. 1.1 is that of Mead and Spitzer [20], 
who experimentally evaluated the value of for various 
semiconductors and enunciated their famous 2/3 E_ law for 

o 

the harrier height in n-type semiconductor of Group IV and 
Group III-V compound semiconductors. In a recent publication, 
Mead et al. [21], have shown that the covalent compound 
semiconductors do not follow the 2/3 law and that ’ ^ ' in 
these cases is directly dependent oil the 'anion' electro- 
negativity of the semiconducting material. Since we are 
mainly interested in Group IV semiconductors, this point 
would not be pursued any further. 

In the Bardeen’s model discussed above the surface 
states have been taken to be localised in space. In 1964, 
Heine [22] proposed a model in which the surface states are 
distributed in space near the metal-semiconductor interface. 
These surface states arise due to the tunneling of metal- 
electron into the forbidden bandgap of the semiconductor. 
Though an empirical formulation of Heine's model was used by 
some authors to study the field dependence of the barrier 
height [23], a detailed study of Heine's model was first 
carried out in 1973 by Pellegrini [24]. The details of 
Pellegrini’s calculations are given in Chapter IV. The 
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horizontal line in Fig. 1.1, connecting Bardeen’s model and 
Heine's model shows that both these models are derived under 
the same premise, that it is the presence of surface states 
which causes the weak dependence of $ B on ^ 

Another model to explain the weak dependence of <J>-g on 
4, ^ has been proposed by Inkson [25]. The basis of this model 
is the reduction in the bandgap of the semiconductor near the 
metal-semiconductor interface. The reduction in bandgap of 
the semiconductor at the surface arises due to the opposite 
effect of the image forces on the conduction and valence 
bands. According to Inkson, the bandgap actually disappears 
at the interface, so that the semiconductor behaves like a 
semi -metal., and the barrier height is determined by the 
energy level at which the top of the valence band and the 
bottom of the conductance band coalesce [2]. 

This brings us to the end of the description of the 
flow diagram given in Fig. 1.1 excepting the small dotted 
block, named the present work, which concerns itself with 
certain ambiguities, to be shortly described, which have 
been observed in metal-semiconductor junctions. These 
ambiguities have been discovered by us in the experimental 
work done by some workers [26] to distinguish between the 
TB, TFE or FE mode of current transport in their fabricated 
diodes. In order to describe these ambiguities and methods 
(as developed in the present work) to resolve them, it becomes 
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necessary to describe the nomenclature and methodology 
developed hy Saxena to distinguish between TE, TEE and PE 
modes of current transport and those developed by levine [19] 
and Crowell [27] to describe the ' T * effect. Hence, in the • 
following two sections we present brief descriptions of these 
works and then elaborate on the ambiguities and methods of 
resolving them. 


1 . 3 Th e Experimental Method to Identify the Mechanism of 
Current Transport in a Fabricated Schottky Barrier 


diode [“14 


The method proposed by Saxena [14] is based on the 
following observations; 

(i) While the TE theory predicts a forward I-V rela- 
tionship for V > -2™ of the form [12], 

I = I a exp (qV/kT), (1.1) 

most of the practical diodes follow the relation, 

I = I g exp (qV/nkT) (1.2) 

where 'n' is called the ideality factor. In eqn. (1.1) and 
eqn. (1.2) I s represents the reverse saturation current, 
given by 

3= qj> -n 

I s = SA T 2 exp(- £?p) (1.3) 

(ii) In some cases the ideality factor 'n' is temp- 
erature independent whereas in some other cases it is 
temperature dependent. Moreover, in the latter cases it is 
commonly found that 
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T 

n = (1 + ■jr - ) (1*4) 

where T Q is independent of temperature [28], 

and (iii) It was shown by Padovani and Stratton [6] 
that for a forward bias (V£ , the current by both TFE 

SL 

and PE theories can be represented as 


where 


I - I g exp (V/V-p) 

V I = Y oo coth (V 00 AT) for TFE 
= V Q0 for PE 


(1-5) 

(1*6) 


V oo = (1-7) 

Prom the above observations, it is seen that the I-V 
relationship under forward bias can in general be written in 
the form, 


I = I g exp(V/Y T ) (1.8) 

where V^ is the inverse of the slope of the In I-V plot. The 
method suggested by Saxena [14] consists of obtaining I-V 
characteristics in the forward region at different tempera- 
tures (T) and plotting V^ versus T, V^, measured at a constant 
current. The typical nature of the resulting plot is shown 
in Pig. 1.2, where the curves X, II and III represent the 
cases where the IE is the mechanism of transport with unity 
ideality factor, with ideality factor greater than unity 
and with I Q effect respectively. Curves IV and V correspond 
to TPE and PE respectively. 
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k T/ q. ( m V ) ► 


Fig. 1.2 Plot of V T Ve rsus ( kT/q, ) to deter minG the 
mechanism of current transport inSchottky 
barriers ( see Text ) . 
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Before concluding, the following points regarding 
Saxena’s work are worth, noticing, 

(i) The experimental method suggested by Saxena [14] 

i 

considers only thermionic emission and not diffusion as the 
transport for emission over the barrier. This was later 
confirmed by the observations of Crowell and Beguwala [20], 
and Rhoderick [31], that Schottky diodes made from fairly 
high mobility semiconductors like Germanium, Silicon and 
Gallium Arsenide etc., should confirm to thermionic theory 
for emission over the barrier rather than to diffusion theory, 
and (ii) the. validity of the experimental methods suggested 
by Saxena [14] should only be tested in very small area 
diodes - this is mainly because of the technological limita- 
tion in getting a large area semiconducting crystal of uniform 
properties. 

1.4 T 0 Effect 

In the previous section on the experimental method to 
identify the mechanism of current transport in a Schottky 
barrier diode, we have already mentioned that many metal 
semiconductor Junctions show a temperature dependent ideality 
factor. Moreover, in some cases this temperature dependence 
has been found to obey the relation 

(1 + T 0 /T) 


n a 


(1.9) 
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These diodes are said to exhibit T Q anomaly, mainly because 
the physical mechanisms responsible for this behaviour were 
not identified in the beginning. It was with the work of 
Levine [19] in 1971, that the first successful model to 
explain this I Q effect was available. Later in 1977, 

Crowell [27] has emphasized that the choice of parameters 
in the model of Levine is not unique, rather several other 
choices are possible. To start with, we would consider only 
the model of Levine which assumes an exponential energy 
distribution, characterized by a characteristic energy E , 
of the interface states. 

Note, however, that the original work of Levine had some 
physical discrepancies, regarding the vanishing of charge 
in the semiconductor and consequently in the surface state, at 
the flat band condition, as highlighted by Crowell, It is 
interesting to note that this physical discrepancy is 
identical to missing a constant of integration in Levine’s 
work. Realization of this fact has allowed us to view the 
models of Levine and Crowell in a novel perspective which is 
presented below. 

1-4.1 E valuation of the surface s t ate charge 

Pig. 1.3 shows the energy band diagram for a Schottky 
barrier made on an n-type semiconductor under (a) forward 
bias, (b) zero bias and (c) reverse bias, with an exponential 




conditions. The amount of charge in each case is supposed to be 
proportional to the dark area D§]» 
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surface state distribution incorporated in it. The negative 
charges will be located in the acceptor like surface states 
which lie above the neutral level [ 19 ] and lie below the 
Fermi -level. The negative charges in the surface state can 
be derived under the following assumptions; 

(i) The surface states distribution in the forbidden 
bandgap of the semiconductor is of the form, 

I) SS (B) = D sq exp (JE|/B 0 ) (1.10) 

where E is the energy measured from the neutral level and is 
positive in the direction towards the Fermi-level, B Q is the 
characteristic energy and D g0 is the density of the surface 
states at the neutral level, and 

(ii) the occupation of surface states follows Fermi 
statistics with zero degree approximation of abrupt transition 
of occupation probability from unity to zero at the metal 
Fermi level. This fixes the integration limits for the 
evaluation of the surface state charge to be 0 and <(> 

where ^ is the neutral level at the surface and ^ is the 
barrier height (see Fig. 1,3). We can write the surface charge 
Q gs as, 

q(<j> -40 

J exp ( I eI / E q ) dE (1.11) 

where |E| •» E, since we are interested only in the positive 
values of E. Integrating eqn. (1.11) , we obtain. 
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Q s s ■= «so E o [exp(4(* H -4, B )/B 0 )-l] 

= Q f [exp(q4 s -* B )/E 0 )-l] (1.12) 

Here- is the interface charge that neutralises, the surface 
state charge at flat hand. We note that the expression for 
Q given in eqn. (1.12) is the same as that used by Crowell 
to modify Levine’s formulation and which marks the starting 
point for Crowell's model [27]. 

1.4.2 Ev aluation of T p 

It is well known that the ideality factor (n) of a 
Schottky harrier diode is given hy the relation [32], 

/ 3 <f> b v "1 . . 

n = (1 - ) (1.13) 

The evaluation of n requires the knowledge of' the voltage 
dependence of the harrier height -g. Prom eqn. (1.12) one 
can write the harrier height ^ as, 

.j, E Q~ ~ t Q_p 

<f> B = « - ~ In (1.14) 

To obtain the voltage dependence of 4 one requires the 

voltage dependence of Q„ 0 . The voltage dependence of CD 

ss ss 

and hence 4 B can be obtained from the charge neutrality at 
the surface as follows. Invoking charge neutrality at the 
surface, one has, 

Q ss “ Q sc (1.15) 
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where Q is the total space charge in the semiconductor 
SC 

and the charge on the metal is taken to he zero. Using 
eqn. (1.14) and (1.15) we have, 



E o Ac + A 
_o ln 


( 1 . 16 ) 


where (1 is given by [51], 

o C/ 


Q sc = f2<le s N d (4, B - 7 - * f "-'TEc7577 
f = r ln = r b 


(1.17 a) 
(1.17b) 


Eqn. (1.16) along with eqns. ( 1.17a) and (1.17b) gives us 
the harrier height as a function of voltage. Differentia' 
ting eqn. (1.16) with respect to V we have, 


As 
a v 


!o 

2 q 


[(♦ 


B 


V - 


klv-l-j / 

q ' J 11 


av 


- 1) 


SC' 
(1.18) 


where f(Q sc ) = 1/(1 + Q f /Q sc ) • 

Erom eqn. (1,18) one obtains the value of (l - 3 ^g/ 3 V) 
required for the evaluation of the ideality factor as, 


^ “ aV 6 ^ = 1 A + 2q ^ B _V “^f " q ^ ’ f ^ 1 * 1 ^ 

o C 

Using eqns. (1.9), (1.13) and (1.19) one obtains 

^r* = (4>-g-V*-4>£- ~p) * f (Q sc ) (1.20) 

* 2E o 

For barrier height <P ^ < <P , when Q gc » Qf » 
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eqn. (1.20) reduces to the following expressions obtained by 
Levine [19]: 



(1.21a) 


= P (In (J/A X T 2 ) - ( b+l) ) (1.21b) 

Q 

where £ -V= ~ln(J/A*l' 2 ) 

Under constant current conditions, ignoring the slowly varying 
In I term one observes from eqn. (1,21b) that T Q is indepen- 
dent of temperature and is proportional to E . 

It is important to note here that the inconsistency 
in the evaluation of surface state charge in Levine's model 
has arisen due to the use of eqn. (1.21a) for T , rather than 
the more general expression given in eqn. (1.20) . 

So far we have looked into the model of Levine and of 
Crowell to explain the T Q effect and the experimental methods 
of Saxena to identify the mechanism of current transport in 
Schott ky barrier diodes. In the following section we present 
a brief outline of the work done in the present thesis. 

1.5 Ihe Bresent Work 

In recent years, due to their applications in photo- 
voltaic [33] and power devices [34], it has become necessary 
to fabricate large area Schott ky barrier diodes. To identify 
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the mechanism of current transport in these cases also the 
method suggested by Saxena [14] has been utilized by some 
workers [26] leading to an ambiguity in this identification 
[35]- The ambiguity is that the experimental plot of V^, 
versus f has shown the mechanism of current transport to be 
field emission for Schottky barrier diodes, made of semicon- 
ductor substrates of low doping [26], when the operation is 
at a high temperature. It seems obvious that in diodes of 
low doping, field emission should not be a mechanism of 
current transport at high temperatures [35]- One of the 
purposes of the present work is to demonstrate that even 
though the mechanism of the current transport in a large' 
area Schottky diode may be TB, the fluctuation of parameters 
can create a situation so that the experimental plot of V«, 
versus T would show BE to be the mechanism of current 
transport. In our study of current transport in large area 
Schottky barrier diodes, we have considered a diode with the 
mechanism of current transport defined by TE with I effect 
and with a fluctuation in E . The choice of the fluctuating 
parameter would be clear from the following discussions. 

1.5.1 Current transport in large area Schottky barrier 

diodes' T 

The choice of fluctuating parameter in a large area 
Schottky barrier diode was based on two important observa- 
tions regarding diodes showing T effect and are given below: 



20 


(i) The first observation is due to Paciovani [32j, who 
noted that the different Au-n-GaAs Schottky barrier diodes 
made on the same wafer showed different values T Q . We have 
also noted a similar trend in different Au~nSi Schottky 
barrier diodes made on the same wafer. 

(ii) One observes from the discussions in the previous 

section that T is a function of E , <J> K and D _ and that the 
O 0 so 

fluctuations in T q could be due to the fluctuations in any 
of the three parameters. To identify the fluctuating para- 
meter, the values of B., <t> and 1) are evaluated using the 
method given in Chapter III, for each diode. It is found 
that while E 0 is different for diodes with different 1' 0 , the 
values of <J> * and D 0 . are nearly the same for all the diodes. 

Guided by these observations, we have considered that 
the value of E Q fluctuates for a large area diode. It is 
further assumed that the variation of areas, AS^ of patches 
with values of E Q between E 0 and E Qi + AE oi has a Gaussian 
distribution [36], i.e. 

AS, 1 (E 0 -E oi ) 2 

= ~rr exp [- a a — 3 aE 0.22) 

5 2 a 2 01 

where is the mean value of E oi ' s and a is the standard 

deviation and S is the total area of the diode. The total 

current flowing in the diode can then be written as, 
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SJ 


= ,;r (<SS i / d E oi )A*I 2 exp(-j^)[eip(g -1] dE oi 


E, 




5 

dE 


dE 


OX 


01 




oil 


(1.23) 


where = the harrier height of the ith path and dS^/dE^ 


is obtained from eqn. (1.22) in the limits AE Qi + 0 as 


as. 


S 


dE - 
01 


V 2na‘ 


(B -E D .) 2 

exp [_ -^-£3-^] 

2o 2 


(1.24) 


To evaluate the value of current using eqns. (1.23), (1.24), 
(1.12) and (1.14), we need the values of V ♦ > C so E oil 
and E oi2 . The area S taken to be 1 cm 2 and the method of 
choosing the values of the other parameters is discussed 
below. 

Keeping in mind that the choice of E Q , 4> * and D &0 
should be close to the realistic values, a preliminary 
experimental investigation has been carried out to choose 
the values of these parameters. The experimental investiga- 
tion consisted of preparing a number of homogeneous Au-nSi 
diodes on a single wafer and evaluating E , 0* and h g0 in 
each diode. The diodes made in the same wafer showed a 
scatter in T Q (measured at the same current density in all 
the diodes) and hence in the corresponding values of E 0 . The 
average of these E 0 's has been chosen to be 1 Q . The limits 
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of integration and ® 0 i2 are ’ taken ^e va ^ ues 

E Q i where the integrand falls to less than 10 % of the maximum 
value. Regarding the evaluation of <j> s and D go , we observed 
that there is no method available in literature to. evaluate 
these parameters. We have developed a method to determine $ * 
and D g0 from the temperature dependence of the barrier 
height at zero bias (see Chapter III). Prom the values of 

and I) _ obtained for different diodes we observed that 
80 

these values were nearly the same for all diodes and are 
0.78 Eg and 1.48 x lO^/cm 2 eV respectively. Here Eg 
represents the band gap of the semiconductor. The experi- 
mentally measured values of $ K and D g0 are then used in our 
theoretical calculations. Using these values and eqns. (1.25) 
and (1.24), the forward I-V characteristic has been obtained 
at different temperatures and the corresponding Ym has been 
calculated from the relation. 


Y, 


I 


at each temperature. 


3 In I _i 

( ) 

3 Y 


(1.25) 


A plot of Vj, vs. kT/q is given in Pig. 1.4 for E 0 =s»075b'V 
with a as a parameter. The following feature of the plot 
should be noted. 

Por a given S Q when a = 0, one gets a plot of vs. 
kl’/q which is typical of T Q effect, i.e. a straight line 
parallel to unity ideality factor line but not passing through 
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the origin. As the value of cr becomes finite a flat region 
appears which is. typical of field emission in the conventional 
interpretation. Thus although the mechanism of current 
transport in the diode is thermionic emission, the fluctuation 
of parameters have brought in the characteristics of field 
emission, if attention is restricted to a limited range of 
temperature. This suggests that caution must be exercised 
while inferring the mechanisms of current transport in 
Schottky diodes from vs. T plots. The analytical detail 
of the calculations of 7^ vs T plot for large area diodes 
showing fluctuations in T Q , is given in Chapter II, where 
we also present a calculation of the effect of fluctuation 
in interfacial layer (oxide layer) thickness in devices with 
a thin interfacial oxide layer. 

An experimental verification of the theoretical calcu- 
lations carried out for large area diodes with fluctuations 
in T is presented in Chapter III, where we also present 
a method to characterise the surface state distribution 
leading to T' 0 effect, and a numerical calculation to establish 
the function of a guard ring in a Schottky barrier diode. 

In Chapter IV we present an analysis to establish the 
i' 0 effect as an effect of metal-electron tunneling into the ' 
forbidden bandgap of the semiconductor and establish an 
equivalence with Levine's model for f Q effect. 
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In Chapter V we present a novel method to determine 
the doping density of the semiconductor in Schottky barrier 
diodes showing I effect. A few concluding remarks are made 
in Chapter VI . 



CHAPTER 2 * 


T HEORETICAL STUDIES OCT THE EFFECT OF FLUCTUATION ON THE 
~ IDENTIFICATION OF THE MECHANISM OF CURRENT TRANSPORT" 


It has been mentioned in Chapter I that the manifesta- 
tion of the mechanism of current transport in the' nature of 
k!T 

the versus ~ plot in a nonhomo gene ous Schottky barrier 
diode is different from that in a homogeneous Schottky 
barrier diode. It should be noted, however, that to demons- 
trate the differences, we have considered only the fluctua- 
tions in 'E ' . It is clear that there would be many other 

w 

parameters which also can fluctuate in a nonhomo geneo us diode 
and affect the identification of the mechanism of current 
transport. The purpose of the study in this chapter is to 
systematically search for the parameters which can fluctuate 
and affect the identification of the mechanism of current 
transport. To do this, first of all, we consider the current- 
voltage equation for a Schottky barrier diode, under a forward 
bias V, (with Y £ *~) [l]s 

I = SA*T 2 exp (-X*6) exp (-ggr^O exp (jjjg?) (2.1) 

Here the term exp(— X^6) represents the transmission coef- 
ficient of the interfacial oxide layer of thickness 6 and 
of mean barrier height X , and other terms are defined in the 
list of nomenclatures. 
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Our aim now is to locate, from a perusal of eqn. (2.1), 
the parameters which may fluctuate and affect the mechanism 
of current transport in nonhomogeneous Schott ky barrier 
diodes. Before proceeding further, let us point out the 
difference between the calculations of current -voltage 
characteristic in a nonhomogeneous diode and a homogeneous 
diode. In the case of a homogeneous diode the calculation is 
one dimensional, whereas for a nonhomogeneous diode it is 
basically two dimensional. The two dimensional nature of 
the calculation in a nonhomogeneous diode arises primarily 
because of the fluctuation of material and geometrical 
parameters transverse to the metal-semiconductor interface. 
These fluctuations may also lead to the flow of a transverse 
current in the nonhomogeneous diodes. To simplify the 
situation, we have incorporated the fluctuations, but have 
ignored the flow of transverse current while building the 
model here. It is worth mentioning here that an experiment 
has been designed (see Chapter III) to check the validity of 
the assumption mentioned above and has shown that under 
practical conditions the neglect of the transverse current 
does not introduce noticeable errors. Hence, the nonhomo- 
geneous diode has in effect been treated as composed of 
several homogeneous diodes, but with varying values of X , 

<f>gQ, n and 6 in eqn. (2.1), in parallel. Furthermore, the 
system has been taken to be isothermal and A* has been 
as sinned to be constant. 
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To look into the possible nature of variation and to 
identify a parameter which, would uniquely define the fluc- 
tuations in x 9 duo* n an<i ^ i " t: is necessar y to consider the 
various barrier height models. From a description of the 
development in the study of metal-semiconductor contacts, 
given in Chapter I, we see that the barrier height models 
can be classified as (i) Work function model (ii) localised 
surface state models or (iii) .Distributed surface state 
models. A study of the effect of fluctuations in diodes 
governed by the work function model has already been carried 
out by Johnson et al. [2], hence is not discussed here. 

The effect of fluctuations in localised and distributed 
surface state models is treated in the present work, the 
former in this chapter and the latter in Chapter IV. 

In the next two sections we present the study of fluc- 
tuations in MIS model (sec 2 . 1 ) and exponential surface 
state model of Levine (sec 2 . 2 ). 

2 . 1 MIS Model 

The study of the effect of fluctuation of parameters 
on the determination of the mechanism of current transport 
in a MIS diode is divided into two parts. The first part 
deals with the identification of a parameter which uniquely 
defines ^ gQ, n and X . The second part deals with the study 
of the effect of fluctuation of this parameter on the 
versus *=■ plot. 
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2-1-1 I dentification of the fluctuating parameter 

For an MIS diode, using the methodology outlined by 
Cowley and Sze [3], and the identification of two kinds of 
surface states by Card and Ehoderick [1] the expressions for 
the barrier height (<J> B0 ) at zero bias and the ideality factor 
(n) can be written as (see Appendix B) 



where ^ = Neutral level at the surface measured from the 

valence band (=0.27 E ). 

°2 “ Vfei + 4 2 «(B as )] 

°1 = 2 « s V 2 /'i S 

B ss = B sa + B sb 

E S a = Density of surface states communicating with 
the metal 

■^sb 88 f en ®ity of surface states communicating with the 
semiconductor.’ 

and D 1 ^ ** Effective density of surface states communicating 
with the semiconductor. 
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Prom equations (2.3) and (2.4) we see that both 4> and n 
are functions of _U ga , P Qb and 6. It is seen from literature 
[l,4j that 6 is the only independent variable and that 1) 
s^d "^sb are ^oth functions of 6. Regarding the mean barrier 
height x of the interfacial layer we see that [1,4] this is 
also a function of the oxide thickness 6. Prom the above 
observations we note that 6 is the only independent variable 
and hence the fluctuating parameter in the MIS model. 


In table 2.1 we present the 6 dependence' of D gg# D gb 
and for 0 < 6 ^ 30 £, evaluated from the data of Card [ 4 ] 
and Card and Rhoderick [1].* Note that we have considered 
the 6 dependence of JD gg and D gb rather than the more funda- 
mental quantities N ga and I) gb because fo.\. the calculation of 4» 

and n we do not need the individual values of h and I) , 

sa s d 

but only those of D gg and IV b . 


BO 


1'able 2.1. Functional relationship between X, D and D' 
and the oxide thickness 5. ss 80 


Parameter 


Function of oxide thickness *6' 


X), 


ss 


D eb 


X* = a6 2 +bj a = 7.10 2Cj eV’^/cm 2 , 


K C6 + dj C 


1.11 x 10 7 eV^ 

= -2.04 x 10 20 /eV cm 3 
d = 6.68 x 10 l3 /eV cm 2 
D gb * e6 3 + ff e = 1.46 x 10' >2 /eV cm^ 

f = 1.28 x 10 12 /eV cm 2 


The value of X) ss at different oxide thicknesses has been 
evaluated by using eqn. (2.3) and a plot of <|> B0 ' versus 6 

given In Ref. r ' n 
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Having identified the oxide thickness (6) as the fluc- 
tuating parameter and defined the relationship of the other 
relevant parameter with 6, let us now look into the effect 
of fluctuations in 6 on the determination of the mechanism 
of current transport in a MIS diode. Before proceeding 
further it should be mentioned that our study is restricted 
to oxide thickness 6 •$ 30 1. 


2.1.2 Eff ect of fluctuations in * 5 1 


The aim in this section is to evaluate the I-V charac- 
teristic of a nonhomogeneous diode, with a fluctuation in 

the oxide thickness, at different temperatures and hence 

kT 

obtain the V (|1 versus ~ plot. To calculate the I-V charac- 

J, H 

t eristic we need to know: (i) the area of the individual 
diodes which make up the large area diode and (ii) the 
temperature dependence of the various parameters, viz. 

Xo» X > £ > Dlv and E„. let us look into these one 

by one. 


In this study we assume that the sub-area A S i of 
patches with oxide thickness between 6^ and dH-Ad.^ has a 
gaussian distribution, and is given by, 


a a i 
8 




( s-d, ) 2 

exp ( * — ) A 6- 

2 a 2 


( 2 . 5 ) 


where 6 is the mean oxide thickness, o is the standard 
deviation and S is the total area of the diode and is taken 
to be 1 cm 2 . 
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i'he parameters jj r , g(1 ^sb and x are all considered to 
be temperature invariant and the temperature dependence of 


* $ m » a 

aid 

m 

X 

are 

given by 

the equations 


yw 

= 

EjO) 

g 

- af 

7(p+i) 


(2.6) 

a 

S3 

7.2 x 

10 “4 

e7/ 0 Kj 

p = llo8°K 



» 

^ mo 

aT 



(2.7) 

^rao 

S3 

4.698 

eV, 

a = 1. 3 

x 10 ~ 4 eV/°K 


X a (‘J-0 

= 

x so + 

bT 



(2.8) 


X so = 4.02 eV } b = 10~ 4 eV/°K. 


The total current I 
analogous to eq. (1.23) as 

6. 


s/V 


s 


f&Ta' 


exp(-x 


flowing in the diode can be written 

9 



exp(~ 


^4 Bois 
"'*kT 1 


exp( Sr } 


exp ( U LzJll ) d6 (2.9) 

2a 2 

where 6^ and are the values of 5 at which the integrand 
drops down to less than 10% of the maximum value. Using 
eqn. (2.3) to (2.9) and eqn. (1.24), the 1-7 characteristic 
and Vrj, under constant current condition is evaluated at 
different temperatures. In Pig. (2.1) we present a plot of 

kl’ ip 

V r j, versus for a fixed value of & for different values of 
a. The following feature of Pig. 2.1 is to be noted. 



6- 

15 A 

0 

cr 

Q.OA 

0 

— , — * cr* 

2 A 

O 

a~ 

3 A 


0*f— 

0 

1 1 
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1 l L 1 1 L 
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ill ( mV) — 1 ► 

Fig. 2.1 

Plot of Vp 

kT 

versus q~ with cr ns a parameter for 


6=15 A- The 

*| ^ . 

donor density is N D ~2-5x1Q /cm 
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The plot of Vj, versus ~ is a straight line passing 
through the origin with a slope greater than unity. The 
slope of the straight line increases with an increase in cr. 
Thus we see that the fluctuations in 6 does not alter the 

wn 

qualitative nature of the 7^ versus ~ plot. 

Prom a preliminary study of the effect of fluctuations 

wn 

on the nature of the 7^ versus plot in the two 

localised surface state models we see that fluctuation in 

kT 

brings about a qualitative change in the 7^ versus ~ 

plot whereas fluctuations in oxide thickness does not cause 

kT 

any qualitative change in the nature of the 7 rj , versus ~ 

plot. Since our interest is in the study of nonhomogeneous 

diodes which do not have a unique relationship between the 

kT 

mechanism of current transport and the versus ~ plots, 

in the next section we shall look into the effect of fluc- 
tuations of E Q in greater detail. 

2.2 Effect of Fluctuations in E Q 

A brief account of the effect of fluctuations in E Q 

has already been presented in Chapter I. In this section 

kT 

we present further details of the 7^ versus — plot 

and present a semi-analytical calculation to identify the 

temperature range where 7m remains independent of tempera- 

* kT 

ture. In Pig. 2.2 w© present 7 rj? versus — plots for 
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Fig. 2. 2(a) Plot of Vj versus kT/q, with cr as a 
parameter for £><070 eV- The donor 
density is 2-5 x 1 0 l5 c m”" 3 






Fig.2.2(b)Plot of V T versus ” with tr os a parameter ^for 
E„ = 0-065eV. The donor density is N 0 -2-5x1o'% 
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different E Q with a as a parameter. Apart from the trend 
explained in Chapter I, the following trends are also worth 
noting. 

(i) the range of temperature over which the value of 
V,£ remains to he independent of temperature reduces with 
the reduction in a for the same value of E . We also 
see that for the same value of cr, the mid point of the 
temperature range over which V t( is constant increases 
with a decrease in E . 

(ii) For a given value of E , there is a threshold 

value of a, below which the V rn versus plot is still 

i Q. 

a straight line not passing through the origin. It should 
be noted however that in the presence of fluctuations in B , 
the slope of the Vj versus — - plot remains to be less 
than unity, as compared to the unity slope in the homoge- 
nous case. 

kT 

Having looked into the nature of the V^, versus 
plot for devices with fluctuations in E 0 , the next logical 
step is to look for a relationship between V^, ' (the value, 
of Ya, when independent of temperature) and T (where 
'ij, ^ I 1 ^ 'i‘ u is the temperature range over which V,^ remains 
to be independent of temperature) and E q and a, to 
identify ohe nature of fluctuation in a nonhomo gene o us diode. 
In the next section we present a semi -analytical calculation 
to evaluate V,^, !T^ and as a function of E Q and a. 



38 


2.2.1 Calculation of T, and. I 

. 


Prom Chapter I, we see that the total current in a non 
homogeneous diode with a Gaussian area distribution in E Q 
is given by, 

h 0 i 2 ■ 2 <{> 

I = SJ = J Sjl exp( ~ j&~ ) exptg) a^. (2.11) 

B oil 


where a[ 


dS. 


oi f2 % 0 


s “( g 0 ..E 0i ) 2 

— exp u± 


2 o' 


,2 


( 2 . 12 ) 


Differentia-ling In I with respect to V to obtain V l£ , we 
have, 


^ = rf -hr + 111 B /° 12 s i “oiJ 

011 (2.13) 


kl 


E oi2 3<i> w 

J S! exp( - -^) iS 01 


"oil 

"-g-~ 

oi2 ~q$ p . 

/ si (-Tar 1 ) dE oi 

^oil 


(2.14) 


In eqn. (2.14), we can obtain 4> B± and (3<f> B /&v) as a func- 
tion of E oi from eqn. (1.14) and eqn. (1.18) of Chapter I, 
as, 


E - 

exp(~^r~) = exp(^||~) (^- c -^ + 1 ) W ~ 


, Bi T oi 

^ TT S TFT- 5^) 


fi 


(2.15) 


Substituting eqn. (2.15) in eqn. (2.14) we get. 
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E 


3 Ini 
3V 


SL . f i 
k'l ^ 


1 

T 


E 


oi2 Q 

S s: (i + 5 


E 


01 

«ci\Sr~ 


oil 


fi 


T . 

01 

(in 
V 01 7 


dE 


E 


E 


oi 2 V- /Tf-rp 

/ s. (1 -I- r—--)"oi //kl - dE, 


oil 




oi 


oi 


■} ( 2 . 16 ) 


fe [i 


where 


p(t) 


T 


E, 


E, 


-012 , klT T . 

, ^ s “ (1 + C Tirffer 


E 


oi2 


0 V 
111 ’ 


; s' (i + 

“oil Qfi 


clE 


01 


(2.17) 


(2.18) 


Our interest now is in finding out the temperature 

range over which Vm remains a constant. This means the 

1 3V T 3 (1/Vj 

evaluation of the temperature range over which — »=■ or — • =— 

3 i 3-L 

is equal to zero. So diff erentiating eq. (2.18) with 
respect to temperature and equating It to zero we have, 

3(1/V 


'■ 3T 


0 


3 

3T 


-Sf ^ ~ 


i.e. [(1 - p(D) ~ + 1-il L . gI ] 


T 1 3 T 


0 


1 * O « 1 


|i(T) 




3 T 


(2.19a) 

(2.19b) 

(2.20a) 


Taking p(T) « 1, we have (l - p(T)) - 1 and the equation 

(2.20a) reduces to 

1 


3joiTl 
3 T 


T 


0 


(2.20b) 
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It is obvious from eqn. (1,14) , which, gives as 

the solution of a transcedental equation and the expression 

for ^sci 5 ^ e( l rL ‘ (1*15).) » that a completely analytical 

solution of eqn. (2.20b) is not feasible. To overcome this 

we have explored to find out a mathematical fit to the 

dependence of (1 + -#■— ) E oi/ kT on 3 Prom a plot of 

, Q sci 4 fi 01 

(1 + q—~~) versus E . given in Pig. (2.3) for a 

**f i 

constant current at different temperatures, we note that 
Q . E . /kf 

(1 + q ) ‘ can be represented by a three piece linear 

ui 

approximation; 

Q an . B^./kT 

(1 + ^7^ = exp C a j E oi + V 1 ')] 3 


( 2 . 21 ) 


with exp(b 3 (T) - b 2 (T)) = AT + 33 
and exp(b^( T) - b 2 (T)) * CT + D 


( 2 . 22 ) 


Using eqn (2.21) and (2.22) in eqn. (2.20a) we obtain, 


z N 1 (AT+B) x x + x 2 + ( CT+h) x 3 
p(T) = rjj • "( aT+B) ' x 4 + x 5 > ( CT-i-hT'i^ 


(2.25) 


where x. 


( E - E .) d T . 

exp . Ta+i^TfT exp( a j E oi> iE oi 


! exp ( 


■k 


2cr‘ 


(B ~E . ) 2 
v 0 01/ 

2 o 2 


1 ,< j ^ 3 
1 ^ k ^ 3 

exp ( aj E 01 )dB oi 

l(k,<3 

4 j £ 6 



1 + Qsci^ E o'7 kT 
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Fig. 2.3(b) Plot of^' ^~/1 Tt versus Eoi at different 
temperatures for 1^*07 5 eV andcr*-025eV . 
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where 


= 

[Soil- E oill] 


H 

5S 

[Soill- B oi21 

and 

h 


[ B oia- s oi2 ] 


are the three intervals over which the three piece linear 

approximation is made. 

Putting, q = Ax-j_ + Cx^ 

p s= Bx, + x 5 + 1 jx~ 

p (2.24) 

s = + CXg 

r = Bx^ + x^ + B'Xg 

and differentiating eqn. (2.23) With respect to temperature and 
using eqn. (2,20b), we have, 

T 3 + T 2 + 1 - pr = 0 (2.25) 

s 2 s 2 

i.e. l l3 + u*T 2 + v'T - w* =0 

with u' = — , v' = r 2 -2ps/s 2 and w*- = pr 

s 2 

Normalising the temperature 1 with respect to 100°3£, we have 
the normalised -form, 

l' 4 3 + uf' 2 + vT* - w = 0. (2.26) 

where T’ = T/100°K. 
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A perusal of eqn. (2.26) with the values of u,v and w 
obtained from our numerical calculation shows the presence 
of a root close to the origin, which can be evaluated using 
the expression, [6] , (Appendix C) , 


, t s= -(v 2 - uw) w/(v 3 -uvw+w 2 ) (2.27) 

It is, however, seen that the temperature corresponding to 
this root falls close to the transition region between 
thermionic emission and thermionic field emission. Based 
on this observation, the cubic equation was reduced to a 
quadratic equation, after eliminating the root 't'. The 
solutions T^ and T of the quadratic equation is obtained 
from, 


T-, + I = u-t 
1 u 

and T 1 T u = w/t . 


(2.28) 


In Big. 2.4 we present a plot of u, t and w as a func- 
tion of a for three different values of E . For the cases 

Jr rp 

where a flatness in versus — is seen, the values of T^ 
and T u obtained from eqn. (2.27) and eqn. (2.28) was found 

to be well matched with the values obtained from the 

kT kT 

versus — — plot. For the cases where in vs ~ did not 

show a flat region, the roots and T u obtained from eqp.. 

(2.28) were found to be complex. 
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Fig*2.4(b) Plot of u ,w and t versus cr"w ith Eo~0*070eV- 



Fig*2.4(c) Plotof u,w and t versus cr it h Ecr 0*06 5 cV. 
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It is important to note that in the temperature range 
T-^ T ^ I the expression in eqn. (2.27) becomes negative, 
which indicates (from eqn. ( 2 . 19 a)) that 3 (l/V^/gT is 
positive. Shis means that in this temperature range the 
value of (l/Vfjr,) would increase with temperature, i.e. 7^, 
would decrease with temperature. It is seen, however, that 
the rate of decrease in 7m is so small in this region that 7^ 
can be regarded to be a constant. In the temperature range 
t < T < and T > T , the value of 7^ increases monotoni- 
cally with temperature. 

Prom the calculations presented in this section we see 

that there is no closed form relationship for the evaluation 

of 7j, T-j. i * 1 terms of E Q and a. This means that the 

value of I Q and a for a nonhomo gene o us diode can be obtained 

led} 

only by template matching of the 7^, versus plot. The 
steps to be followed for template matching are as follows. 

(i) Compare the experimentally obtained 7^ versus 
~ plot with the theoretically obtained plots for different 

M. 

B 0 with a as a parameter. If the experimental plot matches 
with one of the theoretical plots, then the values of E Q 
and a are directly obtained as that of the theoretical plot 
which matches with the experimental one. In case the 
experimental plot does not match with any of the theoretical 
plots, it would then fall between the plots for two different 
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values of a for a certain value of' E . This would, give us 
the value of E Q and the value of a can he obtained by the 
method given in the next step. 

(ii) Having obtained the value of using the rela- 
tion + 1’ = (u-t) x 100, the value of a can be obtained 
from the corresponding plots of u and t as a function of a. 

As an illustration let us consider the plot of 

Trip 

versus ™ given in Fig. (2.5) for a Au-nSi diode made in our 
laboratory. Comparing this plot with the various theoretical 
plots, we observe that it falls between the plots for 
a = 0.035 eV and a = 0.025 eV with E Q = 0.075 eV. This- 
gives us the value of E Q to be 0.075 eV. From the experi- 
mental plot of V r j\ versus kl/q we obtain the value of T^+T 
to be 728°K. Using the relationship T-+T. = 100 x (u-t), 

-L Li 

we get (u-t) = 7.28. From Fig. (2.5b) wherein a plot of u, 
w and t versus a is presented for E Q = 0.075 eV, we obtain 
ct = 0.032 eV for u-t = 7.28. Thus, we seo that for the 
nonhomogeneous diode whose versus kl/q plot is given in 
Fig. (2.5) the values of E Q and a are 0.075 e¥ and 0.032 eV 
respectively. 



In this chapter, a study of the effect of fluctuations 
on the V,ji versus f plot for a nonhomogeneous Schottky barrier 
diode has been carried out to identify the cause of ambiguity 
in the determination of the mechanism of current transport. 






Fig. 2. 5 A Plot of Vf versus kT/q,f or a large area diode 
with N D = 2-5 xlO 15 / c m 3 
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To cLo this the approach taken here is to look at the situa- 
tion as one-dimensional with the fluctuation of parameter in 
the other dimension. This basically means that the non- 
homo geneo us diode is considered to be composed of small 
homogeneous diodes, each identified by a particular value of 
the fluctuating parameter. 

To identify the fluctuating parameters two kinds of 
diodes are considered viz. (i) diodes with a thin inter- 
vening oxide layer of thickness '6' between the metal and 
the semiconductor (Tunnel MIS diodes [18]), and (ii) diodes 
showing T q effect. The fluctuating parameter in the MIS 
diode is identified to be the oxide thickness ’6'. In the 
diodes showing I effect, . the fluctuating parameter has been 
identified to be E 0 » the characteristic energy defining the 
exponential distribution of the surface states. 

It is seen that the ambiguity that arises in the ' 

determination of the current transport from the versus 
kl 1 

— =• plot can be explained through the fluctuations in E . 

It has also been shown that the ambiguity has arisen due to 

the measurements having been carried out in a restricted 

range of temperature. One way to resolve this ambiguity is 

kT 

to check the V T versus — behaviour over a larger range of 
temperatures, extended preferably in the lower temperature 
range. The importance of the extension of the temperature 

* 
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range of measurement in the lower temperatures will be clear 
from the following discussions. 

Extension of the temperature range of measurements 

only to the higher temperatures would not resolve the 

ambiguity, as a monotonic increase after a flat region in 
IcT 

the V,ji versus plot also represents the case of thermionic 
field emission being the mechanism of current transport [7]. 
On the other hand, a monotonic increase in with tempera- 
ture at lower temperatures followed by a plateau, possibly 
occurs only in the presence of fluctuation of parameters in 
a nonhomogeneous diode. 

Having presented the theoretical calculations in the 
effect of fluctuations in T 0 , we shall go over to the experi- 
mental verification of the same in the next chapter. 



CHAPTER 3 


EXPERIMENTAL stupy of the effect of fluctuations in e q 


In the previous two chapters , we have postulated and 
theoretically demonstrated' that in the case of Schottky 
barrier diodes the temperature dependence of the ideality 
factor is affected by the fluctuation in some parameters. 

In this chapter an experimental verification of this fact is 
given. This has been done on metal (Au)-n Silicon Schottky 
barrier diodes which have been fabricated in our laboratory 
for this purpose. 

The basic aims of this investigation are the following: 

(i) It has been observed by Padovani [1] that several 
small area Au-nG-aAs diodes made on the same wafer exhibited 
different values of T . It should he recalled that this 
observation had provided the starting point for the theore- 
tical model given in the previous chapters. Thus, before 
this theoretical model could be applied to any system, it 
must be experimentally checked that T was different on 
several small area diodes made on the same wafer. The first 
aim in this chapter is to experimentally investigate whether 
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such, fluctuations take place in case of metal -n silicon 
diodes. If the answer is in the affirmative then one 
proceeds to the step given below. 

(ii) It is obvious that the theoretical analysis of a 
device in which the parameters are different at different 
points of area would require a two-dimensional analysis 
because transverse currents can flow. Only in the case that 
the role of the transverse current is not dominant that a 
simpler theoretical analysis using one -dimensional equations 
but with fluctuations superimposed can be done. It should 
be recalled that we have used the latter approach in the 
theoretical model developed in Chapters I and II. Thus, we 
have inherently assumed that the role of the transverse 
current is not important in explaining the temperature 
variation of the ideality factor. Our interest now is in 
checking whether this assumption holds true in the case of 
metal-n Silicon diodes. We have designed a simple experi- 
mental procedure to check this point, rhe basic idea is as 
follows. Prepare a large number of isolated small, area 
diodes and measure the I 0 in each. If the 'i' Q ’ s are different 
then connect a few of the diodes in parallel and measure the 
temperature variation of ideality factor. If this plot shows 
features similar to the temperature variation in an identical 
system but with a large area then obviously the role of the 
transverse current is not dominant and the assumption made 
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earlier would hold good. Otherwise a two dimensional theore- 
tical model should he developed. 

The details of the experiment and the results are given 
later (sec. 3*5 and 3.4) after describing the fabrication of 
bhe diode (sec. 3*1 and 3.2) . Since che leakage current of 
a Schottky harrier diode hinders the interpretation of the 
experimental results, we start the discussion with a study 
of the guard rings which suppress the leakage current. 

3 * 1 Guard - r ing Strictures 

A major problem in the case of Schottky barrier diodes 
made on low-doped semiconductors arises due to the presence 
of large surface leakage currents. The surface leakage 
currents are comparable to the normal diode currents, which 
would flow under reverse bias and small forward bi as . Thus, 
it becomes obvious that for any meaningful interpretation of 
current voltage relation in these cases, the surface currents 
must be suppressed. The method commonly used to suppress 
or reduce the leakage current is to employ guard ring struc- 
tures [1]. In published literature three kinds of guard 
ring structures are used, and are shown in Fig. 3.1. 

The first type shown in Fig. 3.1a uses an extension of 
the metal film over the passivating oxide layer. This is 
referred to as the MOS guard ring (IiOSGR) structure. The 
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(a) lb) (c) 

BACK METAL CONTACT 
SPACE CHARGE REGION 

Fig.3.1 The three kinds of guard ring structures used in 

literature. ( a) Mos guard ring ( b) p-n junction guard 
ring (c) Auxiliary Schottky barrier guard ring. 




Fig. 3.2 A cross-sectional elevation of a Schottky barrier diod 
with an auxiliary Schottky barrier guard ring 
structure . The* rectangle abed marks the boundar 
of the region over which the Poisson's equation fs 
solved. 
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MOSGR acts as a field plate to keep the underlying surface 
in depletion when the diode is either reverse biased or 
slightly f oi ward biased, thus reducing the excess currents 
at the corners [2]. The second type is shown in Pig. 3.1b 
and uses a diffused p-n junction guard-ring. Similar to the 
MOSGR here also the guard ring keeps the underlying struc- 
ture in dejiLetion, thus reducing the currents at the 
corners [3]. The third kind of guard ring structure is 
• s * 10wn 3.1c and uses an auxiliary Schottky barrier 

diode surrounding the main diode at a small distance, 
usually of uhe order of 1 mil [1], In this structure the 
two dimensional nature of the potential distribution is 
exploited to obtain the guard ring action by biasing the 
auxilliary diode in such a manner as to form a space charge 
channel between the main diode and the auxiliary diode, 
thus reducing xhe field at the edges of the main diode [4]. 

I)ue to the ease of fabrication and the simplicity of 
the measurement set up for an unencapsulated device, we have 
used MOS guard ring structure. It is important to note that 
the yield using this type of a guard ring is poor. To 
overcome this, we have fabricated a number of devices and 
have chosen those devices which satisfied the following 
conditions: 

(i) The current density under a reverse bias of 10 V is less 
than 5 mA/cm 2 at room- temperature. 



(ii) The diode does not show time -dependent I-V characteristi 
i.e. the meter readings do not show any appreciable 
fluctuations . 

This gave an yield of 25 .percent for roughly one hundred 
diodes fabricated. In cases where the yield of the device 
assumes importance , one will have to go either for a diffused 
p-n junction guard ring structure or an auxiliary Schott ky 
barrier guard ring structure. 

Normally the thickness of the passivating oxide layer 
used in the MOSGR structure is of the order 2000-4000 2. 

In the present case we have used an oxide thickness of about 
200 2 - 300 2 and the reason for the reduced oxide thickness 
will be clear from the following discussion. 

It has been shown by Yu and Snow [5] and Kano et al. [6] 
that when the underlying surface in the MOSGR is accumulated 
or inverted the guard ring action of the MOSGR fails and it 
contributes towards an increase in the surface leakage 
current, rather than suppressing it. In our case, we are 
thus interested in increasing the accumulation voltage and 
also see that in the range of forward bias used the under- 
lying surface remains to be depleted. Based on the observa- 
tion made by liar [7], that the flatband voltage in metal- 
Si 02 ~n-Si system decreases with an increase in the oxide 
thickness, the MOSGR guard ring is used with an oxide thick- 
ness much less than the thickness normally used in literature 
[5,6]. Note that the accumulation voltage obtained in our 
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laboratory for an MOS diode (Au-nSi system) with, an oxide 
thickness of 200-300 2 at room temperature is found to be 
in the range of 0.7V to 1.0V, which goes far enough into 
the forward bias region not to affect the forward I-V 
characteristics. 

Before proceeding further, let us digress a little 
and briefly look into the design guidelines for an auxiliary 
Schottky barrier diode guard ring. 

3.1.1 Design guidelines for an auxiliary Schottky barrier 
diode guard' ring structure 

A perusal of the literature on the auxiliary Schottky 
barrier diode guard ring structure shows that, though it has 
been extensively used by several authors [1,4,8], the spacing 
of the auxiliary diode and the main diode is chosen very 
arbitrarily. In this section we attempt to give a design 
guideline for the spacing between the main diode and the 
auxiliary diode. 

I‘he choice of the spacing between the main diode and 
the auxiliary diode is governed by the effectiveness of the 
auxiliary diode in modifying the field at the edges of the 
main diode, i’his means that to choose the proper spacing 
between the two diodes, a knowledge of the field profile for 
different spacing is necessary, 'i'he field profile is 
obtained by solving the two-dimensional Poisson's equation 




with, the proper boundary conditions. 

* . 1 

The Poisson's equation, considering the mobile car- 
riers (majority carriers) and the charge due to ionised donors 
can be written as [9], 

liz + 3_V 41f li ri — qv\ 

jj? if = sy t 1 - (3.1) 

Normalising the potential with respect to the band bending 
given by (-T + Og- $ f ) where V is the applied potential, 
the distances with respect to the depletion width, and the 
charge density with respect to (qN D ), we have, 

8 2 #5 3 2 a 

rp + - etp(-t s f )] (3.2) 

where * s = q(~y + - * £ )/ac. 

Note that Nig. 3.2 gives the cross-section of a Schottky 
barrier diode with an auxiliary Schottky barrier diode.- 
The boundaries for the solution of the Poisson's equation 
marked by the rectangle abed. 

Considering the half width of the auxiliary diode to 
be greater than thrice the 'depletion width* under all bias 
of interest (see Appendix C) , we can write the boundary 
conditions (refer Pig. 3.2)). 

(i) Along ad 

»(0, y) =1.0 0 < y< A ( 3 . 3a ) 

and from the continuity of potential and displacement vector 
we have, 



4> (0 + ,Y) = <J> (0 ,y) a < y $ A' 

3 4" c 4 

il 4>(0 ,y) ~ c 3 -~ 4>(0~,y) =0 A < y A' (3.3b) 

where n = e /e . 

b O 

(ii) Along dc 

Prom symmetry conditions 

fy 4>(x, A') =0 (3.4) 

(iii) Along be 

<f> (B, y) =0 0 4 y ^ A' (3.5) 

(iv) Along ab 

Prom symmetry conditions, 

0) =0 Ofx N < B ( 3.6) 

The details of the solution of eqn. (3.2) with the 
boundary conditions given in eqn. (3.3a) to eqn. (3.6) is 
given in Appendix D, and b-;re we present only the results 
of the computations. It should be emphasized that these 
results are in normalized form, hence are of general value. 
The normalization factors of all the relevant parameters 
are given in Table 3. 1 . 

3?ig. 3*3 gives the constant electric field magnitude 
profile under reverse bias in the absence of a guard ring 



Table 3*1. Normalisation factors used in the solution of two dimensional 
Poisson's equation. 






structure, whereas Fig. 3.4 gives a plot of the constant 
electric field magnitude profile for three different 
spacings at a constant bias (zero bias with^g = .8 volts). 
The following features of the figures should be noted; 

(1) In the absence of a guard ring structure (Fig. 
3.3) 'the field at the edge of the metal gets to be as high 
as 8 to 10 times the field at the interior. 

(2) In the presence of a guard ring (Fig. 3.4), we 
see that as the spacing is reduced the edge field decreases 
due to an enhanced interaction between the main diode and 
the auxiliary diode. This in turn reduces the edge current, 
thus improving the I-V characteristics of the main diode. 

JDe normalising the distances we see that, for a base resis- 
tivity of 30 Kbhm-cm, a spacing of 0,3 mm [4] would show a 
guard ring effect from zero bias onwards, whereas for a 
base resistivity of 10 Kohm-cm, it would require a spacing 
of about 0.15 mm. Based on the above observations we see 
that lower the resistivity of the semiconductor base, 
lower would be the spacing required to achieve the guard 
ring effect. 

To sumnarise, one observes from the plots given in 
Figs. (3.3) and (3.4), that for an effective guard ring 
action, one requires the spacing at the maximum forward bias 




Bg.3.4(a) Constant field profile near the edge of a Schottky 

* barrier diode with a guard ring at a distance of 3»2W . 




Fig. 3.4(b) Constant field profile as the ed 
diode w‘ith a guard ring at a d 







beyond which the leakage current is negligible, to be less 
than 2.5 times the depletion width at that bias. 

Having looked into the choice of guard ring structure 
for the diodes under study, we now go over to the details 
of fabrication and the experimental setup in the next 
section. 

3.2 Details of Fabrication and the Experimental Setup 

In the first part of this section we present the 
details of diode fabrication in our laboratory and in the 
second part a description of the experimental setup is given. 

3-2.1 Details of fabrication 

In Pig, 3.5 we present a flow diagram of the fabrica- 
tion process for Au-n Silicon Schottky barrier diodes. The 
following points about the fabrication is to be noted: 

(i) All the chemicals used in the cleaning of the silicon 
wafer are of electronic grade, (ii) Initially an oxide of 
1500 & thick is grown, which is later reduced to about 
200-300 % after cutting windows in the oxide layer. This 
is necessary because better results are obtainable in photo- 
lithography with thicker oxides. After the completion of 
the fabrication process, the device is taken (without 
encapsulation) to the measurement table. 



Fig.3.5 A Flow diagram showing the diode fabrication proce 



71 


3.2.2 lie tails of the Measurement setup 

The details of the measurement setup are shown in 
Fig. 3.6; Fig. 3-6a shows the block diagram of the setup 
for the current-voltage (I-V) measurements and Fig. 3* 6b 
shows the setup for the capacitance-voltage (C-V) measure- 
ments. The meters used for measuring voltage and current 
are standardized with respect to a standard source and the 
error bar is found to be about +1$ over most of the range 
of interest in the 1-7 measurements. In the .IV range of 
the digital multimeter the error bar was about +3$. 

The variation of the temperature of the sample from 
80°K onward is obtained by using an arrangement shown in 
Fig. (3.7). Several choices were considered for obtaining 
a variation in the operating temperature of the device 
under test. The major constraint under which the choice 
had to be made was that the samples did not have bonded 
leads, hence a probe had to be used for making electrical 
contact to the front metal contact. The main difficulty in 
providing for the movement of a probe on cooled sample 
ensues from the fact that moisture quickly gets deposited 
on cooled exposed sample. To obviate this difficulty the 
system is designed in such a fashion that a positive pres'- 
sure of dry nitrogen was always maintained near the sample. 

The setup consists of a copper rod of 1 inch diameter 
held inside a copper tube of 1*5 inch diameter. The inner 



2 




Fig.3.6(a)Block diagram of the I-V measurement setup, 



Fig.3.6(b) Block diagram for the C-V measurement setup. 
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copper rod is used for mounting the sample under test and 
carries a few windings of Nichrome heater coil at the top 
to vary the temperature. The copper tube with the copper 
rod is suspended from a bakelite sheet mounted on a slotted 
angle frame. The copper tube is fixed to the bakelite sheet 
at one end by means of a brass flange attached to the copper 
tube. Apart from supporting the copper rod, the tube also 
serves as a jacket around the copper rod, through which a 
continuous flow of dry nitrogen is maintained. The dry 
nitrogen is introduced into the copper tube through a small 
hole in the tube just below the bakelite sheet. 

The portion of the copper tube carrying the copper rod 
which is below the bakelite sheet is introduced into a 
Dewar. The other end is closed with a bakelite lid, after 
loading the device under test. The bakelite lid has a 
central hole to introduce the contact probe. The central 
hole in the lid also functions as an outlet for the stream 
of dry nitrogen. Before filling up the Dewar with liquid 
nitrogen, the whole system is first completely flushed with 
dry nitrogen to remove any trace of air in the system. 

The temperature variation is obtained by heating the 
copper tube to varying degree. This set up gives us a 
controlled temperature from 80° K onwards, correct to +2°K. 
The temperature of the system is measured at the top of the 



rod using an Iron-Constantan thermocouple. The thermo- 
couple was standardized at liquid nitrogen temperature, 
melting ice temperature , room temperature and at boiling 
water temperature, with the reference junction kept in a 
melting ice bath. 

This brings us to the end of the description of the 
fabrication method and the measurement setup. In the next 
section we present a method for characterizing a Schottky 
barrier diode showing T Q effect. 

3 . 3 Characterisation of a Schottky Barrier Diode Showing 
T Q Effect 

Before presenting the result of experimental investiga- 
tion, let us consider the methods for the evaluation of some 
of the parameters of interest in a homogeneous diode exhibiting 
T q effect. These parameters are E Q , the energy characterising 
the assumed exponential distribution of surface states of 
the equivalent Levine’s model, the barrier height at 

zero bias* , the neutral level at the surface measured 
from the bottom of the conduction band, and I) . the surface 
state density at the neutral level. Methods for the deter- 
mination of each of these parameters are given below. 

3.3.1 Determination of E 0 

The value of E Q can be determined both from the forward 
and the reverse I-V data [10]. We present the method of 


determining I) 
° o 

one by one. 


from the reverse and the forward 1-7 data 


(i) Evaluation of from the reverse 1-7 data: 

Under a reverse bias 7 >, ~ — , the reverse current I x 


' q. ’ 


can be written as. 


Ir = A*T 2 expfq ({i-g/kl 1 ) 


(3.8) 


Prom eqn. (1.14), for Q gc » Q^, we have 

♦ „ = ** - f° in (2a£) 

V B v q v Q f ' 

hence we can rewrite eqn. (3.8) as, 

* nA* E /kl -E_/kT 

I R = AT 2 exp (- ||“) (Q sc ) (Q f ) J 


_ E_/kI 

I* (tj o/ 


-E /kT E /kT 


where I* = A*T 2 exp (^^r) (Q f ) ° ( £ ) 


(3.9) 


(3.10a) 


(3.10b) 


(3.11) 


Q sc ' M\ 

— = If— (# B - V " *1 ~ — > 


(3.12a) 



(-7 + a) 


(3.12b) 


where a 




B f 4 


Prom eqns. (3.10b) and (3.12a) we see that a plot of In 
versus In (~7+a) will be a straight line of slope m 2 =(E 0 /2kT) 
and the value of E Q can be obtained from this slope using 
the relation. 
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E 0 = 2kT m 2 (3.13) 

It should he noted here that in our studies we have taken a 
mean value of a = 0.4V at all temperatures. 

(ii) Evaluation of E Q from the forward I-V data: 

A method of evaluating E Q from the forward I-V data 
has been proposed by Levine [10] and is as follows: 

Prom eqn. (1.21a) we obtain the voltage dependence of 
2? 0 of the form, 


fj* ( *3) “ v - *t ~ (3.14) 


1_ _ 2£ / d> A k! 

T o ' 

This equation predicts that a plot of ir” versus V would be 

1 o 

a straight line (taking <J>g to be a slowly varying function 
of voltage) with a slope (m) given by, 


m = 


% 
3 T 


2q 

e 0 Tt+t o ) 


(3.15) 


Note that the slope 'm* has been obtained from eqn. (3.14) 
using the relation, 


3 V ~ T + T (3.16) 

o 

Let the value of T Q and m evaluated at a mean value of ? be 
called and m respectively. 1'hen we can obtain E Q by 
using 5 and in eqn. (3.15). 

We suggest that if a plot of 1/T Q versus lnl is used, 
rather than the 1/T Q versus V plot, then it is easier to 



determine E . The reasons for this can be seen from eqn. 
(l.21b), which, is reproduced below? 

= - |^ [In (l/A*T 2 S) + b + 1] (3-17 a) 

= - |£ [In I + b + 1 - In (a*T 2 S)] (3.17b) 

o 

where I is the current flowing in the diode. Prom eqn. (3.17a) 

1 

we see that the plot of ts— versus In I will be a straight 

o 

line of slope m^ = -2k/E Q . 'therefore, one can obtain, the 
value of E Q from the slope m^ using the relation, 

s 0 = (3.18) 

where k is the Boltzmann's constant. 

The advantage of using the versus lnl plot compared 

1 ° -i 

to the 7 s— versus V plot lies in the fact that while 7 jr~ versus 

0 1 0 
V is approximated to be a straight line, the 7 p— versus In I 

i o 

plot is an exact straight line. 

3.3.2 Determination of <5 

Bo 

The barrier height at zero bias, $ cart be obtained 

from (i) a plot of versus V, (ii) ™ - V plot and 

o c 

(iii) the extended saturation current from the forward I-V 
relationship, let us briefly look into these methods one 
by one . 



(a) Prom 1/1' versus V plot : 

By requiring the plot of versus V to he approximated 

1 o 

hy a straight line with an intercept Y^, it has been shown 
by Levine that the mean value of is less than <f>^ the 
barrier height at the voltage V ± , by an amount [10]. 

Thus, <J>gQ can be obtained from the intercept of the i— 

l o 

versus V plot as, 




Y. 

i 


+ $„ + 


kT 

q 


E o 

2Q 


(3.19) 


(b) Prom the forward I-Y plot : 

It has been shown by Levin [10] that the forward 
current I, at any voltage Y can be written of the form, 


I se exp(qV/k(T+T 0 )) 


( 3 . 20 ) 


where I_„= the intercept of the tangent at a current I 


'se 


on the current axis 
= SA*T 2 exp (~q<J> B 0 A(l , + , l o ) ) 

Prom eqn. (3.21) we see that we can obtain the barrier height 


(3.21) 


4* BO aSj 




BO 


k(T+T ) 

[In (I se /0A S 1 2 )] 


( 3 . 22 ) 


(c) Prom the 1/C 2 -V plot : 

For Schottky barrier diodes showing 1 0 anomaly, it 
has been shown by Levine [19] that the capacitance C 



80 


can be written in the form 

HIM*# 


C 2 


1 2SZ— ( x + _2. + _£./4 x) 


T3 2 

r 

w 

2 ' 


where 


x = (4> 


B 


Y 


q. 


ki\ 
q *' 


(3.23) 


Por x large, i.e. small forward bias and reverse biases we 
have, 


1 


E 


^ir ( hn - V- r- - * „ - S + -a.) (3.24) 

S U 1) 


c 2 ~ e a qK Tl v BO ' T+T ” f q q 


, . 3<!> „ T 

where ^ — 9 j. _ tt _ <J> 4 . tt , . . Q 

wnei e B BO 3 Y v BO V f+T 


(3.25) 


(J 5 " Egq^S 2 Ti+T 0 /T) Ca+V f )(*B0-V q + q ^ ~ V ] 

(3.26) 

where T Q is the average value of T in the range of voltage 

over which the C-V measurement is made. Prom eqn. (3.26) 

we see that barrier height <j>.gQ can be obtained from the 

intercept Y^ of the ~ _ y plot using the relation, 

C 


❖ 


BO 


V, 


T 


(T+f 0 ) 


+ K + 


kT E o 


(3.27) 


3.3.3 .Determination of <J> and 1) 


so 


A perusal of literature shows that there are various 
methods to determine $*, the neutral level at the surface 


measured from the bottom of the conduction band and D the 

ss 
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surface state density for Schottky barrier diodes with, a 
thin interfacial layer [11,12] whereas for devices showing 
T effect, there is no method to determine <J> * and D . It 
should be noted that in the case of devices with an inter- 
facial layer, the energy distribution of the surface state 
density (D gg ) is considered a constant, on the other hand, 
for devices showing effect Levine's model assumes an 
exponential energy distribution of surface states, given by 
the equation (ref. Pig. 1.3) 

E ss (£) = 2 S0 exp [(E-q«fr*)/E 0 ] (3.28) 

Prom eqn. (3.28) we see that for the surface state density 
D gs , we require, apart from d * the values of D g0 and E Q . 

A method to evaluate E Q has been given earlier and here we 
will be presenting a method to determined and D . 

SO 

We see from eqn. (1.14) that the zero bias barrier 
height at any temperature T can be written as 

Q ( 1) 

<!> bo CI) = **(t) - E 0 In [-SS2 + 1 ] ( 3 . 29 ) 

where Q f = q ’D QQ E Q 

= space charge in the semiconductor at zero 

S CO 

bias. 

Since this is the only equation available for the determina- 
tion of ** and D sq , we need to look for one more relation to 
obtain d* and This we obtain from the temperature 

SO 
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dependence of 4 >gQ, under the following assumptions. 

(i) We assume that is given by the relation, 

<f, s (T) = r E g (T) (3.30) 

where E (T) is the hand gap of the semiconductor at a 
e> 

temperature I! and p is a proportionality constant. The 
temperature dependence of the band gap for silicon is 
obtained from the relation [ 11 ], 

E g (T) = E g (0) - aT 2 /( (3+20 (3.31) 

where E_(0) = 1.16 eV 

s 

a = 7.02 x 1CT 4 eV/°K 

and p = 1108 °K. 

(ii) Normally, Q sco > the zero bias space charge in the 
semiconductor, is a temperature dependent quantity, but by 
restricting ourselves to the temperature range where Q 

° SCO 

is a constant, we see that the plot of ^^(f) versus E ( T) 
will be a straight line of slope r . 

Under these assumptions we see that the intercept of 

IS Q 

the 4> B0 ( T ) versus E g (T) plot gives ~~ In ( + 1 ) and the 

slope gives r . Thus 41 36 at any temperature T can be obtained 
from eqn. (3«30). Knowing the barrier height ^ B q(T) at a 
fixed temperature, we can calculate Q sco and hence from 
the intercept 4 > B0 using the relation, 
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Q-p — Q 

1 SCO 


[exp(q^ 0 /E o ) - l] 


-1 


(3.32a) 


Prom eqn. (3.18) and the relation Q„ = q D E we obtain 

1 SO 0 


^ g0 from the equation, 


1) 


so 




(3.32b) 


o ^O “O 

Thus D g0 can be obtained once the value, of <j, B0 is obtained 
experimentally as a function of temperature. This concludes 
our discussion of the methods for determining E Q , <J> B0S <> * 
and ^so for a diode showing T effect. 


3 * 4 Results and Discussions 

In this section we present the In I-V and ~ - V plot 

C 2 

at different temperatures, and other relevant plots required 
for the characterisation of (Au-nSi) Schott by barrier diodes, 
as derived from the In I-V plots. 

The basic trend of the results obtained is best des- 
cribed through the characteristics of four diodes made on 
the same wafer. Three of these diodes (EAu 110, EAu 120 and 
EAu 130) are of 1 mm diameter, and the fourth diode (EAu 410) 
is of 4.2 mm diameter. Pigs. (3.6 to 3.12) give the In I-V 
plot of EAu 110, EAu 120, EAu 130, EAu 410 and CD 123, the 
'diode* obtained by the parallel combination of EAu 110, 

EAu 120 and EAu 130. Pig. (3.13) gives the V^, versus 
plot for EAu 110, EAu 120 and EAu 130 and Pig. (3.14) gives 
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Fig, 3, 14 A plot of Vj versus kT/q for the large ar< 
diode EAu 410 and combined diode CD 123 



kT 

'Che I/™ versus — — plot for EAu 410 ana Cl 123. The following 
**■ Si 

features of the plots in Pig. (3.13) and. (3.14) are worth 
noting. 

(1) From a plot of versus —■ given in Pig. (3.13) for 
the three diodes of 1 mm diameter, we see that in all the 
three diodes the current transport is by thermionic emis- 
sion with T effect. The values of T for EAu 110, EAu 120, 

0 o 

and EAu 130 are respectively 27°K, 38°K and 48°K. Note that 
T q was measured at a constant current of 1.0 mA. 

(2) The plot of Vm versus --for EAu 410 given in Pig. 3.14, 

-*• Si 

shows a constant value for V ^ for a temperature range 270°K $ 
T £ 356°K and for temperatures T < 210°K and T > 356°K 
Vfj, increases monotonically with temperature. Hence we see 
that though the small area diodes surrounding a large area 
diodes exhibited thermionic emission as current transport, 
a measurement of the large area diode restricted to the 
temperature range 250°K T ^ 380°K, would lead one to an 
erroneous conclusion that the mechanism of current transport 
in this case is field emission. 

and (3) Prom the V^ ( versus ~ plot of the combined diode 

CD123, we see that though the, individual diodes showed 

thermionic emission as the current transport with T Q effect, 
kT 

the versus plot of the composite diode goes through 
a temperature range (275°K $ T ^ 325°K) wherein v^, versus 
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T'C ^ P o (~) 

— - flattens out. The value of fox 1 T < 275 n and. T > 325 K 
is seen to increase monotonically with temperature. 

kT 

Thus we see that the versus — plot of the parallel 
combination of three isolated small area diodes with dif- 
ferent T 0 's shows features similar to that of an identical 
system but with a larger area. This means that the trans- 
verse current does not play a dominant role in explaining 
the temperature variation of the ideality factor in a large 
area diode and hence the assumption made earlier while 
developing the theoretical model holds good. 

Having verified the basic assumptions and the results 
of our theoretical model, we now look into the characterisa- 
tion of a diode showing T Q effect and identify the parameter 
whose fluctuation results in the fluctuation in T 0 . The 
following plots which are necessary for the characterisation 
of the three small area diodes are given in Pigs. (3.15) to . 
(3.20) : 

(i) Plots of 1/T q versus V are given in Figs. (3.15) to 
(3.17). 

(ii) Plots of ^£q( t ) versus E^T) for all the three diodes 
are given in Pig. (3.18). 

(iii) Plots of In I R versus In ( -V -I- 0.4) at room temperature 
for all the three diodes are given in Pig. (3.19). 

(iv) A plot of 1 /T q versus In I for all the three diodes 
is given in Pig. (3.20). 



356 K 



V ( volts) — *■ 

Fig.3.15 Plot of versus V at different temperatures for EAulfO 
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Fig. 3 .18 A plot of <f>g 0 (T)vBrsus Eg(T)for the three small 
area diodes EAu llO, EAu120 f E Au 130 . 



Fig. 3.19 A plot of In I r versus ln(-V+ 0*4 ) at room temperature for the 
three small area diodes E Au 110, EAu 120, EAu!30* 
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O, 0 

In table 3.2 we present the values of E o s ^BO ^ at different 

temperatures), 4> s and J) for each of these diodes. There 

so 

we have also included the results of two other sets of 
diodes made on different wafers. 

The following points are worth noting from our results. 

(1) The harrier height determined using all the methods 
mentioned in sec. 3.3, match well with each other. 

(2) The value of B , obtained from the forward 1-7 data 

o 

and the reverse 1-7 data are very nearly equal, and 

(3) Bor a given wafer, different diodes with different T Q , 
show a variation only in E Q , while <{> and D remain 
constant. This means that most probably E Q is the 
parameter, whose fluctuation causes the fluctuation 
in I 

So far we have been cqnsidering the nature of the I-V 
characteristics and the various characteristics derived from 
them. We now more over to study the effect of T q effect, 
and E q fluctuations on the C-7 characteristics. 

Before we present the experimental plot of - 7 for 

C 

the diodes made by us, let us look into the basic difference 
between the — 7 plot for a tunnel MIS diode [18] and for 
a diode showing T anomaly. It should be noted that the 
basic motivation behind this comparitive study is to look 
for an effect of T Q anomaly on the - 7 plot as opposed to 
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a thin, interfacial layer ‘between the metal and the semicon- 
ductor and to use the difference, if any, to identify 

1 

diodes showing T Q anomaly from the - Y plot. 

C 

(i) Diode showing T effect : 


We see from eqn. (3.19) that ig- is given by 

c 


si 

c 2 


1 


E„ 




kT, 

— H 

L 

(3.26) 


¥e see from eqn. (1.21a) that for measurements made under 
constant voltage conditions (T /T) is a constant [10] and 
hence the average value (5 Q /T) will also be a constant, 
evaluated over the same voltage range at all temperatures. 
Thus, we can rewrite equation (3.13) of the form, 

F = if -f + r ) - T] (3 - 33) 

From eqn. (3.21) we see that the effect of T q effect on the 
- V plot is two fold and are : 

c 

(1) The slope of the ~ - V plot is scaled down by a factor 

0 . . 

np= (1 + T /T) and 

E 

(2) i'he intercept is given bjr, = n^V^ + ~) , where 
is the .intercept in the ideal case. 

(ii) Diodes with a thin interfacial layer ('funnel MIS diodes) 
We see from the basic definition that, 


rt uu „ j ___ 

§ - - “at 2 = tv [V2qe s N D ^ B - V - <J> f - TEETaTTl (3.34) 
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where Q = total space charge in the semiconductor 

O v 

C = depletion capacitance 
and V is negative for reverse "bias. 


Carrying out the differentiation eqn. (3.34) reduces to, 


C 

1 


* £ s 


2(4 b - V- 4 f - ) '9V 


84 

(; 


B 


1) 


(3.35) 


3 ($) *n ~| 

Using the relationship n = (1 , eqn. (3.35) is re- 

written as, 


q e N,, 
H s U 



S 


2(*B ■ 

- T ”*f 

kl\ 



1 

l 

2S" 2 

(A •tr 


kT 

* V • f 

c 2 

n 2 


(«> B - V 

- 4 f ■ 


Using 

the 

relation, 4^ 

= ^B0 + 

34 p 

W 

V, i 


(3.36) 

(3.37) 
V, eqn. (3.34) reduces 


to 


— 2 


1 2S r / i . kl\ 

0 2 - qSgNjj ' n '[ n ^B0 ” *f “ 


(3.38) 


From eqn. (3.38) one observes that the effect of an inter- 
facial layer is to scale up both the slope and the intercept 
of a — - V plot, as compared to an ideal diode. 

' (y ■ ■ ■ ■ 

Comparing eqn. ( 3. ?:6) and eqn. (3.28) we see that the 
only difference between the ^ - V plot of diodes with an 

c 

interfacial layer and diode showing T q effect is that in 
the former case the slope of the ~r - V plot underestimates 

n <• 
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the doping density, while in the latter it is overestimated. 

The basic nature of the — - V plot in both the cases is the 

C 2 

same and the differentiation between the two models from 

the - V plot is not possible. In Pigs. 3.21 to Pig. 3.23 

we present a plot of versus V at different temperatures 

for EAu 110, EAu 120 and EAu 130. In Pig. 3.24 and 3.25 

1 

we present a plot of — — versus V at different temperature 

c 

for EAu 410 and CD 123. The following features of these 
plots are of interest. 

(1) The plot of —y versus V is a straight line for small 
area diodes with a constant slope at different temperature 
as predicted by eqn. (3.26), and 

(2) The nature of the ^ versus V plot for a large area 

Lf 

diode and the diode formed by the combination of the three 
small area diodes do not show any considerable variation from 
that of a small area diode. 

3.5 S ummary 

In the previous two chapters we have developed a theo- 
retical model to show that in the case of Schott ky barrier 
diodes the temperature dependence of ideality factor is 
affected by the fluctuation of some parameters. The theore- 
tical model is derived based on the following assumptions. 

(1) There is a scatter in the values of T Q for the different 
homogeneous diodes made on the same wafer. 







Fig. 3*2 2 t/c 2 versus V plot at different temperatures 
for the diode £ Au 120 . 
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Fig- 3-2 5 c2 versus V plot at different temperatures 
for the combined diode CD 123. 
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(2) The theoretical analysis of a device in which the 
parameters are different at different points of area would 
require a two-dimensional analysis, as transverse current 
could flow. We have reduced the two-dimensional analysis 
to a one-dimensional analysis superimposed with the fluc- 
tuation in the other dimension, with the assumption that the 
transverse current does not play a dominant role in the 
nature of temperature dependence of the ideality factor. 

In this chapter we have set about to experimentally 
verify the above mentioned assumptions. To do this following 
simple experiment has been carried out. A large number of 
isolated small area Au-nSi diodes are fabricated and the 
f Q in each of these diode is measured. It is found that 
different diodes made on the same wafer have different 
values of T Q . A few of these diodes are connected in 
parallel and the temperature variation of the ideality factor 
in the combined diode is measured. This plot is found to 
show features similar to the temperature variation in an 
identical system but with a large area, thus verifying the 
assumption that the transverse current does not play a 
dominant role in the temperature dependence of the ideality 
factor in a large area diode. 

The next step in our experimental study was to estab- 
lish a parameter, the fluctuation of which results in the 
fluctuation of T Q . To do this, we require the values of 
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^ sxid u so , the three parameters which characterise the 
exponential surface states responsible for the 1* anomaly. 

A perusal of literature shows that there is no method 
proposed for the determination of the values of <J>* and jj 

so 

in a diode showing T Q effect. Here we have developed a 
method to determine the values of a s and f from the 
temperature dependence of the zero bias barrier height 
The method outlined to determine ^ * and u go is formulated 
under the following assumptions. 

(1) The neutral levels®, measured from the bottom 
of the conduction band bears a constant ratio T to the 
band gap 3^(1) at all temperatures. 

(2) The total depletion charge Q gco of the semicon- 
ductor at zero bias is nearly a constant over the tempera- 
ture range of interest. 

• * ' ■ • ‘ • . / .. ■ . ■/ ' ■ ' \ ' 

/ • " ' ; 

Pfom une values of B Q , and jJ. S0 determined for dif- 
ferent sets of diodes made on the same wafer, we observed 
that for diodes with different T Q 's, the values of B Q 's are 
different while and f go remain constant. This vindicates 

the choice of E Q as the fluctuating parameter in the theore- 
tical model. 

To conclude, we see that our experimental investigations 
justify che assumptions under which the theoretical model 
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lias been developed to study the effect of fluctuations in 
geometrical and material parameters in the temperature 
dependence of the ideality factor in a Schottky barrier 
diode. 



CHAPTER 4 


T Q AS AN EFFECT OF TUMBLING OF METAL ELECTRONS INTO 
THE FORBIDDEN BANDGAP OF THE SEMICONDUCTOR 

Depending upon the physical configuration of the Schottky 
harrier diodes several theoretical models for explaining the 
T 0 effect have been developed [1-3] . Some of these models, 
along with the one that we have developed using Heine's 
model of a Schottky barrier, are the subject matter of this 
'chapter. We start by briefly mentioning the main features of 
each of these models. 

(i) M-N contact with localised surface states : 

A model to explain the T q effect from the temperature 
dependence of the ideality factor in the forward I-V charac- 
teristics was first proposed by Levine [1]. Levine showed 
that for T q effect to take place, the surface states at the 
metal-semiconductor interface should be exponentially distri- 
buted in energy with a characteristic energy centered at 

an energy level <J> *. Thus the expression for the surface state 
charge Q gg that must hold (as shown by Crowell [2]) is: 

Q ss = Q f [exp(-(<j>* -4 > b )/E q ) - 1] (4.1) 

Here (j> B is the barrier height and Q f is a constant charge 
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Starting from (4.1) a relation between E Q and T Q can be 
obtained. It is important to note that according to Levine, 
there is a unique relationship between I and E Q . An 
important conceptual advancement in this area came with the 
work of Crowell, who clearly showed that the relationship 
between 1 Q and E Q was not unique, as a matter of fact, it is 
possible that some other distributions in energy (of surface 
states) may also show T q anomaly provided eqn. (4.2) is satisfied 


al * Q S c 

3 $b 


= constant 


(4.2) 


at all temperatures and bias [2], In eqn. (4.2) Q repre- 
sents the depletion charge in the Schottky barrier diode at 
a given bias. Based on eqn. (4.2) Crowell showed that T 

o 

anomaly can be obtained in the case of parabolic surface 
states also [2]. 


(ii) MIH Contact : 

Based on Levine’s exponential surface state model [1], 
Rhoderick [3] has proposed a model for 2? 0 anomaly in MIR 
contacts. He has modified eqn. (4.1) as follows, 

Q ss + Q m = Q f l> x P ~ ^ jj) /E q - 1] (4.3) 

where Q m is the charge on the metal. 
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(ill) M-P-W contact s 

A model to describe T q anomaly in M-P-N contacts was 
proposed by Crowell [2], who showed a? Q anomaly to be a result 
of the modification in the potential profile of a metal- 
semiconductor diode, due to an inversion of semiconductor 
type at the surface. Pig. (4.1) shows a schematic c har ge 
profile and energy band diagram for a Schottky barrier 
diode with conversion of semiconductor type at the surface 
[ 2 ]. 

(iv) M-N contacts with distributed surface states * 

A perusal of literature on M-H contacts with distri- 
buted surface states [4-6] shows that (to the best of our 
knowledge) no theoretical model for 2? 0 effect has been 
developed for this kind of contacts. In this chapter, an 
attempt has been made in this direction. Por this purpose 
one needs a model of a Schottky barrier diode with distri- 
buted surface states. A model for such diodes has already 
been developed by Heine [4] and detailed calculations using 
this model have been done by Pellegrini [6]. We have 
heavily drawn on the works of Heine and Pellegrini while 
developing our model which is presented in sec. 4. 3. 2. 

Before proceeding further we would like to bring out 
the basic difference between Heine’s model [4] and other 

surface state models [1,7,8]. Phis is done in the following 
section. 
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4 . 1 Generalised Approach to the Surface State Models 

The initial studies on surface state models were car- 
ried out by Bardeen [7] and Heine [4]» Their attempt was to 
explain the experimentally observed weak dependence of 
barrier height on the metal work function in some metal- 
semiconductor diodes. This weak dependence on work function 
could not be reconciled with the direct dependence of 
barrier height on work function as predicted by Mott [9]. 

The basic methodology followed in all the surface state 
models is to solve the Poisson’s equation in the semicon- 
ductor with the condition of charge neutrality providing 
one of the boundary conditions. It should be noted that the 
various surface state models have arisen due to the different 
approaches adopted to evaluate the surface state charge, 
let us now consider the charge neutrality equation and 
identity each of the surface state models. 

For a metal semiconductor contact in the presence of 
an interfacial layer, the charge neutrality condition can be 
written as, 

+ « 8 c + Q ss + = 0 (4.4a) 

where Q ffl = charge on the metal 

Q__ = space charge in the semiconductor 

OW 

Q gg •= surface state charge 
sind Q qx ss charge in the interfacial layer. 
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In this chapter we deal only with the cases where Q ox = 0, where 
the charge neutrality elation simplifies to the relation 


The starting point for Cowley and Sze > s [8], Levine's [1] 
and Heine's [4] model is equation (4.4) (because they all 
have Q ox =0). In the models developed by Cowley and Sze [8] 
and Levine [1], both of which are/onl&deen's surface state 
model [7], the surface state charge (Q ss ) is taken to be 
localised at the semiconductor surface, m Heine's model 
W Q ss 18 distributed in the semiconductor bulk close to 
the semiconductor surface. For the evaluation of Q^, Cowley 
and Sze's model [8] considers a charge, free interfacial 
layer of thickness '6', across which there is a potential 
drop 4 . Choosing a • Gaussian-pill' at the metal-interfacial 
layer interface, assuming that there are no localised 
charges at this interface, one can write the charge in the 
as equal to the flux flowing out from the interfacial 

layer, 

^out Q m = e ox^t) X (4.5) 

where e ox = di-electric constant of the interfacial (oxide) 
layer 

%>x ~ field in the interface layer. 

to opposed to Cowley and Sze's model [s], both Levine' s [1] 
and Heine's [4] model assume a direct contact between the 
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metal and semiconductor resulting in a situation where one 
would have to solve Poisson's equation in the metal with the 
proper boundary conditions to obtain the charge in the metal. 
In the case of Levine's model [1] the calculation of Q m is 
got over with by assuming that Q m = Oj while Pellegrini [5] 
in his calculations based on Heine's model, obtains the 
charge on the metal by assuming that the potential in the 
metal (v ) can be described by the relation, 

v m ( x ) = v 0 exp (-x^) • (4.6) 

where v Q is the potential at the metal -semi conductor inter- 
face, and X is ' tlie Thomas Fermi screening length in the 
metal. 

From Fig. (4.2) wherein we have presented the energy 
band diagram for the various barrier height along with the 
respective potential profile for the Poisson's equation, we 
note the following: 

(1) In the case of localised surface state models, the 
presence of the surface states govern only the boundary 
condi tions*a*» Vhile the charge profile and hence the nature 
of the potential profile remain to be the same as in the work 
function model, and 

(2) In the case of distributed surface state model, the 
surface state charge not only modifies the boundary condi- 
tion, but also the space charge term in the Poisson’s 
equation. 




FJg^.2 A schematic ef Hit various barrier height models with the 
corresponding energy bond diagram and charge profile. 
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Having considered the various surface state models 
and their essential differences, let us next consider the 
models to explain the effect. 


4.2 Existing Models for T Q effect 


The various models advanced to explain T Q effect can 
broadly be divided into three categories, (i) Work function 
model, (ii) Surface state model, (ill) poped interface 
model. let us now look into them one by one. 


4-2- 1 ^|^ ° a r model piOT fltf. contact free 

This model was advanced based on the following obser- 
vations [10], 

(i) In the case of diodes showing I Q effect, . the reverse 
saturation current follows the relation. 


where S 
^ mso 


saV exp {4 mso A(aM-T 0 )] (4.7) 

area of the diode, and 

effective barrier height (independent of 
temperature). 


and (ii) The temperature dependence of the barrier height 
follows the temperature dependence of the work function ($ ) 
Of the metal, the electron affinity (X g ) being assumed 
temperature independent. 
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Prom eqn. (4.8) the barrier height 4 g(T) at temperature I 
is given by, 

* B< T > = * m /' 1 + V T > (4 ’ 8) 

Considering the temperature dependence of the metal work 
function of the form, 

V 1 ’ = 4 mso - aI < 4 -» 

and using the relation 0 „ = $ -X we have, 

±5 in s 

T o * T C^mso /cc:r) " 1 ]" 1 (4.10) 

The work function model fails because of the following 
reasons: 

(i) It does not explain the bias dependence of the 
barrier height observed in most diodes [10]. 

(ii) It does not explain that the value of T q , when 
evaluated under constant current conditions, becomes inde- 
pendent of temperature. 

4-2.2 Surface state models 

The barrier formed between the metal and semiconductor 
in the presence of surface states can be classified into 
strongly pinned barrier or weakly pinned barrier. Strongly 
pinned barriers are a limiting case of the surface state 
model, where there is a large surface state present on the 
semiconductor surface [11], The weakly pinned barriers are 
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the harriers which fall between the two extreme cases of 
strongly pinned barrier, where the barrier height is indepen- 
dent of the metal workfunction and the workfunction model 
where the barrier height is directly proportional to the 
workfunction of the metal. We shall be looking into the I 

anomaly xn the case of weakly pinned barriers and the strongly 
pinned barriers separately. 


Before proceeding further, let us consider the neces- 
sary condition for a temperature independent T Q [2]. Assuming 
that the image force lowering is negligible and that the 
potential energy maximum is defined right at the junction we 
can express the space charge Q sc as 


2N r . i 

« so = (-S (r, - 1) S)» 

s * 


(4.11) 


where 




*( % - v) 
0 


ln 


*c effective density of states in the semi- 


conductor conduction band. 


For an applied forward bias V > we have 


Jin I 


JL 


and 


7 ] 


~ + T 0 ) ~ E0 

= In (A*T 2 /J) - in (W C /W I) ) 


From eqn. (4.11), (4.12a) and (4.12b) we obtain. 


(4.12a) 


(4.12b) 



kf 

o 

q 
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3 <*) _ js . 

— — ~ V ss — ( 3 <J>_. /3 Q ) • (3 Q /3- f] ) — — • 7 — • . ... - 

5 P V B' ^sc ; v ^sc 7 n ' 2 (rPXJ T^lnCl 

s c 

(4.13) 

Thus if T q is measured at constant q it will be temperature 
independent only when (3<Jw3lnQ ) is a constant over the 

-D SC 

bias range associated with the operation at constant q. It 
should be noted that operating at constant q, if the logarith- 
mic variation of temperature in eqn. (4.12b) is negLected, is 
equivalent to operating at constant current. 

Having considered the necessary condition for a temp- 
erature independent T , let us look into the models for T 

U Q 

effect in weakly pinned barriers and strongly pinned barriers 
one by one. 

(1) T Q effect in weakly pinned barriers : 

There are three models to describe T effect in 

o . 

weakly pinned barriers and are (a) Levine's exponential surface 
state model [1], (b) Crowell’s parabolic surface state model 

[2] , and (c) MIS model with exponential surface states [3]. 

Let us now consider them one by one. 

(a) Levine’s exponential surface state model. 

The first model for T Q anomaly in Schottky barrier 
diodes proposed by Levine [1] has already been dealt with 
in sec. 1.4 of Chapter I. Here we bring out only the salient 
features of this model starting from surface state c ha rge of 
the form* 
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Q ss = Q f [exp (q($ £ -$ B )/E 0 ) - 1] (4.14) 

We have from eqn. (4.14) with Q gs = Q gc (considering Q m =0) 
we obtain. 


E. 


♦ B - ♦* ( 


Z t *t 


■) 


where 


” ' t { 

€ = Q ctr /£ i md SL = Q f /e_. 


(4.15) 


Using eqn. (4.13) and eqn. (4.15) we obtain 

1 




2lV>iT * T't+ %) * 


where 34» B /9ln£ = £ • E Q /(g > + ) 


( 4 . 16 ) 

(4.16a) 


Under the condition that for a weakly pinned barrier 'T/t-^ >> 1> 
we can obtain, 

rp" 1 _ 2k_ (n - 1) 
x o ~ qE 0 lT} ±j 

= (* B - v -<fr f - l^)" 1 (4.17) 

Using the expression J = A*T 2 exp (q(V - 4^) /kT) , where ^ 

is the barrier height at the voltage ’ ~ is re-written in 

^o 

terms of the current density J as, 

Ir = P * [In ( J/A S T 2 ) - ( b+l)]” 1 (4.18) 

o^o * 

where b = q ♦ = kT In (Ng/N^) . 

Erom (4.18) we observe that under constant current conditions, 
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due to the very slow variation on InT, T Q can be considered 
temperature independent. 

(b) Crowell’s parabolic surface state model [2] 

Prom eqn. (4.15) the variation in^g'is given by, 

-f + 'tp g -t 

1<V $ B0> = B o ln = B o 111 [1 + fU% ] U - 19) 

where <J>gQ is the barrier height at a field =^ 0 - Taylor 
expanding eqn. (4.19) up to the second term, and neglecting 

£ ~t Q 

the higher order terms (^ — << 1) eqn. (4.19) reduces to 

^•o u 

fd <J> ) E ° ~^°1 i ~^ > 1^ (A on) 

' B “ BCr ~ q [^ o+ ^] - 2 [ -£ +£ f J (4.20) 

Equation (4.20) defines the barrier height -fie Id relation 
in the parabolic surface state model. In Pig. (4.3) we have 
reproduced from Crowell [2] a plot of (4> ^ • - 4> fi ) /g versus 
(S -€ 0 )/(g +-6j) for a comparison of the exponential 
surface state model and parabolic surface state models. 

Pig. (4.3) shows a very good match between the plots 
representing the parabolic and the exponential surface 
state models. 

(c) MES model 

T q effect in MIS model was identified by Rhoderick 
[3] by including the metal charge in Levine's foraulation. 
The charge neutrality in the presence of a finite metal 
charge is given by, 




< t-U)/(t»+%) ► 


Pig.&*3 Comparison of electrical field dependence of the 
harrier height in units of ^o^ac cording to the 
exponential interface state model and two para- 
bolic interface models, one of which matches value, 
slope and curvature af?«* $;and the other which 
features a three point match/ 
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Q ss = Q ss + Q m = %c <4' 21 > 

Following the methodology used by Levine [1], (see sec. 1#4, 
Chapter I), it can be shown that T q anomaly can be seen in 
MIS diode under the condition, 


Q ss = Q m + Q ss = «f [«P(<l(*“ - * B ) /= 0 ) - 1] (4.22) 
with i Q = T 0 Q bo / 9»)/(»Q 8b /»* b ), whereo = (4 B -V) (4.25) 
In the case of an MIS diode Q m is given by the equation. 


Q. 


e ox ^ox 


(4.24) 


m ox '-'ox ox 6 

where is the field in the charge free interfacial oxide 
layer) 6 is the thickness of the interfacial layer) and A 
is the drop across the interfacial layer. From potential 
profile given in Fig. (4.2), we obtain 


hence 


^m 

*1 


e ox ^ m X s “ ^ 


ox 

nr 


(4.25-) 

(4.25a) 


Using (4.24a) into equation (4.24b we obtain, 


T, 


T 


3 Q 3 Q e. 

sc , f ^ss hv 
/ W a “ T“ / 


34 


B 


(4.26) 


30 



1 2b 


(ii) 2? o effect in strongly pinned barriers : 

As compared to a weakly pinned barrier, in the case 
of a strongly pinned barrier, the barrier height is very 
close to <>*, the barrier height at flat band, and hence t f »E. 
Under this condition, eqn. (4.16a) reduces to ' 

3( ft B _ E o 

3ln£ “ “ ~Y^~ (4.27) 

We find that eqn. (4.27) does not satisfy the criterion for 

T 0 effect, hence in a strongly pinned barriers T Q cannot be 
observed. 



Based on the necessary condition for the e sistence 
of T Q anomaly , Crowell [2] has proposed a model for T 
anomaly by considering a doped interface layer between the 
metal and semiconductor. Pig. (4.1) gives a schematic 
charge profile and energy band diagram for a Schotiky barrier 
diode with conversion of semiconductor type at the surface. 

ihe maximum field ^ and the apparent barrier height 
for an M-P-H contact are given by, 




f q( w 

J r 


(4.28) 



(4.29) 
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In equation (4.28) ¥ represents the depletion width in the 
n-type semiconductor. From eqn. (4.28) we obtain for 


N a (x) = ^.(x./x) 2 . 


qN n x 2 . 

= UBl a 


(i 


m 




iV,5 

£ s T 


V (1 


WA 3 ) 2 


(4.30) 


taking ¥/x^ » 1, we have 


e s ” - Q gc (4.31) 

Using eqn. (4.28) to (4.31) we obtain, 

S *B « B D x i r 

jlnQ sc 5l nt e Q ^ ~ (4#32) 

qU D x 2 

= = e q for x m « x. (4.32a) 

From eqn. (4* 32a) we see that T q anomaly can be obtained 
from an M-P(thin) -U contact with the doping in the p-layer 
U A (x) = NqCxj/x) 2 . It should be noted that, whils the doped 
interface model gives a physical origin for E , the model 
is heavily process dependent, as well a£ the metal used for 
the . formation of the Schottky barrier [2]. It should be 
noted -chat this model is restricted to a limited types of 
doping profile in the p— layer. 

4,3 T o aS ^ Effect of tunneling of Metal-Electron into 
■ Forbid d en bandgap of the Semiconductor 

From the various models for T Q effect described in 

the previous section, the following important observations 
are made. 
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(i) The necessary condition for the existence of 3? 0 effect 
in a Schottky harrier diode is that the slope of the 

plot should he a constant at all temperatures and biases. 

* ' 

(ii) Levine [1] showed that to obtain a constant slope of 
lnQ s £*4>B plot, the surface states should he exponentially 
distributed in energy. According to Crowell [3], however, 
this restriction of exponential energy distribution of 
surface states is not the only one that would lead to T 

o 

effect. Por example, Crowell has shown that the potential 
profile of a Metal -P( thin) -N diode is such that one can 
obtain linear lnQ sc -$B relationship for certain doping 
profiles of the P-layer. 

(iii) Heine [4] has considered metal-semiconductor contacts 
with the surface state charges due to the tunneling of metal 
electrons into the forbidden band gap of the semiconductor. 

It is interesting to note that the potential profile in 
this case has a very large correspondence with that of 

M-P( thin) diodes. Hence it seems worthwhile to investigate 
whether in these (Heine's) cases also a T Q effect can be 
obtained. Let us reiterate that I effect can be observed in 
M-P-H diodes for very limited types of profile in the P-layer. 
To search for a similar restriction, if ary, in case of 
Heine's model [4] is the main purpose of the study presented 
in this chapter. The steps to be followed are summarised 
below, 
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(1) Pellegrini [6] has carried out the quantum mechanical 
calculations for Schott ky barriers obeying Heine's model. 

His method can be used (as shown in sec. 4-3.1) to calculate 
the depletion charge Q and the barrier height <f> at dif- 
ferent bias for various temperatures. 

(2) Then one can check for the linearity of the InQ -4-r, 

sc B 

plot at room temperature. This study is repeated at 
various temperatures and the temperature range in which 
lnQsc-^g Plot is linear is noted. 

(3) The forward I-V characteristics is calculated and 

hence the value of various temperatures under constant 

current condition is obtained. Check is then" made whether 
the V rj} versus T plot exhibits a T q effect, in the tempera- 
ture range determined in (2), for the cases where slope of 
lnQ sc -4>2 plot is a constant independent of temperature. 

The details of the calculations are given below. 

4.3,1 Pellegrini's Calculatio ns for Heine's model F61 
X Heine-Pelle gr ini model! 

When a metal is brought into contact with an n-type 
semiconductor a region on the latter to the right of the 
interface (see Pig. 4.2) is depleted of the conduction band 
electron and at the same time its forbidden energy gap 
becomes occupied by electrons as a result of the quantum 
mechanical tunneling from the conduction band of the metal. In 
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an analogous way the electrons of the semiconductor valence 
land may penetrate below the conduction band of the metal 
giving rise to Thomas-Permi charge. It should be noted 
that When the Thomas-Permi charge at the metal is considered 

ay from the metal. This is in addition to the dipole due 
to the charge due to the ionised donors, whose field is 
directed towards the. metal. The presence of these two 
dipoles gives rise to the potential profile shown in 

FXS ' U ’ 2) ' L8t US now oonsilier the salient features of 

the calculations carried out by Pellet-rim- -r* u ■ 

^ -^ ej - -*- e gnni for Heme's model. 

(i) Calculation of the surface state charge 

It has been shown earlier that in the case of Heine's 
distributed surface state model a Pledge of the space 

distribution of the surface state charge is necessary, since 
it contributes to the charge in the Poisson's equation 
Parker et al. [5] have solved the Poisson's equation by 
assuming a surface charge distribution of the form 


p ss = ex P (-x/d) 


and choosing an arbitrary value of * to fit their experi- 
mental data. Pellegrini in his calculations differs from 
ers in that he treats the metal-semiconductor contact 
as a single system and solves (i) Schro dinger* s equation 
(n) Poisson’s equation self consistently for the 


(4.33) 
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evaluation of the space charge due to the tunneling of 
metal electron into the forbidden bandgap of the semicon- 
ductor and the potential distribution respectively. For 
this the potential energy term in the Schrddinger equation 
is assumed to be a constant equal to the potential at the 
metal -semiconductor junction. Using the wavefunction 1 
the space charge due to the surface state is found, which 
is used to solve the Poisson's equation. Thus an improved 
potential energy form is obtained. The process is continued 
till self- consist an cy is obtained and hence the final charge 
and the potential profiles are obtained. Since the investi- 
gations carried out here are of an exploratory nature, we 
will not go into the details of Pellegrini's calculations, 
but only present the results relevant to our study. 

Before proceeding further we would like to point out 
the salient feature of the Heine-Pellegrini models. The 
most important aspect of Pellegrini's calculations is the 
introduction of a phenomenological parameter E^ , the limit 
value of the transverse energy beyond which the electron 
quancum penetration into the forbidden energy gap is zero. 
The genesis of Ej would be clear from the following discus- 
sion. 

Pig. (4.4) shows a plot of E ks -k^ gb for different 
values of E^,, where E^ a is the energy in the semiconductor 




Fig.4.4 Electron energy diagram in the semiconductor 
as a function of kJsb ond Ef 
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and t 2sb ls the propagation vector in the band b of the 
semiconductor in the direction of tunneling. Prom the 
igure, one can see that for a given value of ^ (say E , , 
the corresponding value of J k M . for . - Bf ks 

‘ ^xsb I for Ej > Em is much 
greater than the value of Ik' . for E , 

• xs hi Ior < Em, so that the 

associated value of the wave function becomes negligible 
““ ^ thOSS “ *isb • Here B, represents 

Z T * ln the direotioa tra ~ to *• - 

l t - a Plot of B^ versus 

xsh , ra her than the conventional E^-k p lot . Ihe 

reason for this is as follows. 

qu d ■ J 361161 ' 3117 ’ “ ear the banSed e e (B c or E ) the energy is 

. , ’ d USUally niUe » v analogy with the 

free electron theory as, 


Ej = 

~ h 2 
2 21 

K , 

xsb 


m* 

m xsb 

where k 2 = 

Isb 

k 2 . 
ysb 

+ k 2 
zsb 


+|fi 2 


k 2 

xsb 


m 


xsb 


( 4 . 34 ) 


a a plot of this is given in F±p- (44) j? 

xj.g. For > -r, ,2 

18 Positive, therefore kis re „ °’ 

exists. ■ por energie ■ 1 « solution 

energies in the forbidden gap, E < E < E k* 

“ DeSatlVe ’ ^ 9 ^ on -Wally damped solutions result ’ [12] . 

Having looked into the salient point in the Heine- 

Pellegrini model, let us now look into +* 

■ LO0K lnto the results of 

Pellegrini - s calculations. 



( ii) Evaluation of 


potential profile: 
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2?he charge due to the penetration of metal elec- 
trons are considered to have two components n sc (x) and n fl (x) 
which occupy states in the bandgap near the conduction anl 
the valence band respectively. These „ glven fcy ^ 
relation, 


n sc^ x ^ = N c ex P (“ x A c ) 
n sv (x > = «v exp (-x/ X t ) 


(4.35a) 

(4.35b) 


where H 0> V * 0 and X v are obtained in terms of the known 

parameters of the metal and the semiconductor and the limit 

energy aj. and ^ are the densities of the metal elects 

tunneling close to the conductance and the valence band and 

X ° X v aPe tlle res Pective mean penetration depths 

contribution of n so < x ) and n^x) to the built in potential 
is obtained as 


(N C X c + K X 2 ) 

t s c c V v' 


(4.36) 

incorporating the charges n sc (x) and n sv (x) into the 

Poisson * s equation along with the charge due to the ionised 
donors we obtain, 

3 2 . 


a -Azi j, , , 

1?^ - t (n sc (x > + n sv (x) - h d ) 


(4.37) 


A Closed form solution for eqn. (4.37) is chtained by assumin 
that there exists a length w of a section of the depletion 
region characterised by the inequalities. 
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v c “ V - v i (w)» k r f/q (4.38a) 

W» w » X c , X v (4.38b) 

where v c = ($ m - X g .- <>£,) is the contact potential difference 
v = applied bias 

v(v) = the potential at the point w, 
and ¥ = depletion width. . 

Under this assumption the field at any point within the 
semiconductor can be attained as. 


£(X) = - f s [N c A 0 exp (-x/x 0 ) + I,»v e *P( -a=A v ) 'V W " x) ] 

(4. 33 ) 

Integrating eqn. (4.39) with the boundary condition that 
v(0) = v Q , one obtains the potential profile. 


qU^x 2 

v(x) = V 2— 

s 


(¥ - |) +2^. [l -exp(^x/X )] N\ 


2 

V V 


+ | [l-exp(-x/X c )]lT c x * ( 4 . 40 ) 

s 

* = /If ^ 

v b = ( v o + v q -'I s ) 

The quantities v^, <£(x), and v(x) are all dependent on 

the voltage (V) through the voltage dependence of € and v , 

o o * 

the field and the potential at the surface. These equations 
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along with the condition of charge balance give a transcen- 

al equation for v Q . The solution of this transcendental 
equation leads to the knowledge of v „,-£(*) md r(x) (for 
details see Appendix D) . 

». iii) Evaluation of Current; 


-or the potential profile obtained as above, a general 
expression for the current is obtained as [5,15], 

i = SA*T 2 t exp (-<j> B /k!) [exp (qv/kl) - 1] (4i4x) 

■ ■ ■ . I ' 

where <J> B = barrier height = (<J> ~ X + v ) 


= maximum value of v(x) 

X = effective transparency of the energy barrier 
to the thermionic emission 


q( v c -v+v m ) /kT 

= 1 + / exp(n)/[i + exp(p(|p|j a ]dr, 

a = | p| /p 

In eqn. (4.42) the function F(p) is defined by, 

. x?(p) 

P(Tl) = ^ (2m Xsc )T x / (ti) -l[v m -v(x)-pk!T]^ dp 


( 4 . 42 ) 


( 4 . 43 ) 


where and x 2 are the turning points where (v ffl -v(x)- 

It should be mentioned here that have restricted 

our calculations to the ran^o n -F j 

ange of doping density and applied 



bias, whore T . ls very close to ^nity. ^ 00nfln . ng the 

investigation to the metal-semiccnductor contacts where the 
current transport is predominantly due to thermionic 
emission. This reduces eqn. (4.43) to the simple form. 


SA T erp( -qj> g/ki l ) [ exp( qV/kl’) _ i] 


(4.44) 

This brings us to the end of the presentation of Pellegrini >i 

results and in the next section we present the results of 
our investigation. 

I 

^ ^ ^ Result s and i)j s cus si on 

In this section we present the results of our calcula- 
tions, based on Heine-Pellegrlni • s model, on Au-nSi contacts. 
These calculations have been carried out to expire the 
possibility of finding T Q effect as . result of tunneling 
of metal-electron into the forbidden bandgap of the semi- 
conductor. as already mentioned, to do this we require to 

obtain the lnQ s0 versus * B , and V f versus M plots, and 
check that , * 

(i) the plot of lnQ gc versus 4 B is a straight line, 

(ii) the slope of the InQ versus <f> • 

1 oC ve ^sus v B plot ls a constant 

independent of temperature, 

(ill) in the case where (i) and (ii) are true, ^ ^ versus 

M/q plots (as obtained from the forward I-V charac- 
teristics at different temperatures) exhibit T 
effect and 0 



(iv) in the case where (i) is not true, the 7^ versus M 
plot does not exhibit 2 ? q effect. 

(1) lnQ sc -4 B Plot: 

I'he plots of lnQ sc versus * B , evaluated from eqns. 
(4.35a) to (4.44) for a Au-nSi diode with the semiconductor 
doping Hj, = 2 x 10 16 /cm 3 are given in Hg. (4.5) for dif- 
ferent temperatures with Ej as a parameter. The parameter 
Ei la the limit value of the electron energy in the semi- 
conductor, perpendicular to direction of tunneling, beyond 
which the electron quantum penetration into the semiconductor 
negligible [6], A perusal of Pellegrini's results [6] 

Show that varies from diode to diode for the same metal- 
semiconductor system. Based on this observation we have 

chosen Ej as a variable parameter in our calculations. The 
following features of Pig. (4.5) are of lnterest . 

(i) ‘i'he lnQ gc versus 4 B plot i s a straight line at all 
temperatures for E^> .03 eV. 

(li) i’he Slope of the lnQ gc versus plot for E^,> 0 .03 eV 
is independent of temperature, thus satisfying 
Crowell's [2] criterion for I Q effect. 

Prom these results, one concludes that according to 

Heine- Pellegrini ' 3 model, the T o effect is expected to 

occur for ^ > 0.03 .7. let us now loch into the 7, ve ra us 1 

Plot to investigate the existence of T 0 effect in these 
diodes. 
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Fig.4.5 (a) Plot of InW versus <p B for E-f-0.078 eV- W-Q S c/q,N 0 
lb) Plot of InW versus for Ey «0*065eV. W=Qsc/qN D 



(volts) — ► <Pb (volts) — w- 

(C) (d) 

PiSI. 4.5(c) Plot of InW versus for Et-0 0452«V- W-Qsc/^Nd 
( d)Plot of In W versus^ for Ef.0.028eV- W Qsc/^N 0 
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(2) V rf versus T plot: 


The V T versus % Plot for the diodes, whose lnQ sc versus 

Pl0t 13 giV8n “ +.5. ±3 given in Pig. 4 . 6 . II note 

the following features of these diagrams: 

(i) Por the diodes with E ' > 0.03 eV, the V,, versus I plot 
is a straight line , parallel to the unity ideality factor 
line but not passing through the origin. Thus, we see that 
Schottky barrier diodes in Heine -Pellegrini’ s model [ 6 ] 
exhibit T o effect for certain values of s'. 

Ui) The value of T Q for different diodes, studied here, 
increases with a decrease in s’. 

(iii) Por the diodes with ^ < 0.03 eV (for example, see 

plot for aj = .027 eV in Pig. 4 . 5 (a)) which do not satisfy 

Crowell's condition for T e-f + k« v ■ 

x o 61 foot, the V^j, versus T plot 

can be fitted to the form Vj, = f- (aTtb), where the slope -a- 

is greater than unity. Thus we see that in the case where 

lnQ so-' 1 ’B Plots are not a straight line of constant slope 

independent of temperature the diodes do not show T 

effect, in the conventional sense. 

Similar calculations were carried out for different 
doping densities and the following observations were made: 

There is an upper and a lower limit on the range 
of doping densities, outside of which T Q effect is not 
observed. The upper limit is fixed by the doping density 
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beyond which the thermionic field emission current starts 

to be comparable to the thermionic emission current, and is 

found to he 5 , 10 ls /cm3. me lower limit arises due to 

the reduction in the ideality factor with the reduction 

xn doping density, in our calculations we see that at a 

doping of about ^ = 10 15 /=m\ the ideality factor falls 
down to about 1.01. 

Ihe explanation for the decrease in 'n' with a 
decrease in the doping density i s as follows: 

It is well known that as the doping density ff. 

decreases the depletion charge Q so <„ fly,) als0 deceases. 

Based on this, we see that in eqn. (D.13) (reproduced below 

for convenience) , used to evaluate the potential a t the 

metal semiconductor interface (v ) 

* o' s 


e * v 
m o 


(q^ v X v ~ ql^W) 


(4.45 

the contribution of the term decreases with decrease 

in the doping density Thls means that ^ ^ valu<j Qf 

““ 13 deCreasea ' «» v Q and hence the barrier heig 

approaches the flat band barrier height. Thus, we see that 

as the doping density is reduced the barrier approaches the 
case of a strongly pinned barrier, wherein the barrier heig! 
changes very slowly with the applied bias, leading to an 
ideality factor close to unity [2]. 
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To conclude, we observe that T Q effect can also be 
obtained when the surface states are due to the tunne ling 
of metal-electron into the forbidden band gap of the semi- 

proio.de d suitable value of sj and doping densities 
are chosen. 


Equivalence of Eft and 7 ^ 

llie characteristic energy E Q in the exponentially 
distributed surface state model has been taken to be an 
empirical parameter. No attempt to link it with physical 
causes of surface state distribution has been made in 
Levine work [1], Having established that Heine-ftsllegrini • 
model, under suitable constraints, can lead to T q effect, 
and that the value of T Q is dependent on the value of e' , 
it becomes tempting to attempt a correlation between E ^d 
If this succeeds, one would have succeeded in lining 
the phenomenological parameter E Q with a much more funda- 
mental phenomenological parameter E^. 

The first difficulty one faces in doing this is due to 
the observation that the surface states in Heine's model 
are distributed both in energy and space (near the metal- 
semi conductor interface) whereas in Levine's model these 
“* distributed in energ, but localised in space. The other 
difficulty arises from the observation that while Levine 
considers a aero charge on the metal, Heine's model considers 
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a finite charge on the metal. To see the correlation men- 
tioned. above, we integrate the surface state charge for the 
He me -Pellegrini model, in space and add the charge on the 
metal and check for the relation [ 3 ] , 


Q sc - Q rn + Q ss = Q t [e^pf-qfc *-* B )/y-l] (4.46) 


If the exponential relation of eqn. ( 4 . 46 ) is obtained then 
the correlation of s' and E Q can be easily established. We 
see that for the case of Qj « Q sc and by taking logarithm 
on both sides eqn. ( 4 . 46 ) reduces to: 


ln<2 = lnQ . 
sc ^sci 


c 


B 


where lnQ 


sci 


= In Q f + Jt- « in q + 

O 1 


r V T) 

E~ 


(4.47a) 

(4.47b) 


Hote that use has been made of the relation $ * = f e (T) whic: 
was developed in sec. 3.3.3 of Chapter III. Using the plot 


lnQ sc versus ♦b* ^^sci at different temperatures can be 
obtained, which in turn is plotted as a function of E (T) 

to obtain the values of Q f and r . Once r is Imown J* is 


also obtained. 


In Table 4.1 we have presented the values of E , $ 

Q f 3113 “so (“so = VV for different sj and N^„ Bote 

that B so represents the equivalent surface state density 

at an energy qd* from the bottom of the conduction band. 

From the table we see that, (i) The values of d* and b are 

. • SO 
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Table 4.1, Table showing a correlation between and E Q 

for three different values of the doping density. 


Doping Density 
ffj) ( cm - 3 . ) 

Bj, (eV) 

E 0 (eV) 

(eV) j 

j 

®so 

( cm ~2 eV - l) 


0.078 

0.022 

.8 E g (T) 

4.94 x 10 11 

2.2 x 10 16 

0.065 

0.024 

.8 E g (T) 

4.92 x 10 11 


0.0452 

0.036 

.79EJT) 

& 

4.94 x 10 11 


0.078 

0.031 

.81E g (T) 

9.20 x 10 11 

6.3 x 10 16 

0.065 

0.035 

.8 E g (T) ' 

9.40 x 10 11 

| ■■ ■ . . ■ 

0.0452 

0.043 

. . . i 

.81E g (T) 

9.40 x 10 11 


0.073 

0.038 

.8 E g ( T) 

2.30 x 10 12 

2.2 x 10 17 

0.065 

| 0.043 

•79E g (T) 

2.10 x 10 12 


0.0452 

0.047 

.791JT) 

2.30 x 10 12 




g 
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nearly independent of the value of E^ for a given doping 
density, (ii) The value of E Q increases with a decrease in 


and that 

the 

values of E Q can be fitted to the 

form. 


E 

0 

= a .+ b/E’ 


where 

b 

= 0.00078 eY 2 

(4.48) 


a 

= [1 + logtEp/lO 16 ) ] x 0.014 eY 



Carrying out numerical calculations to evaluate the depletion 
charge Q gc at zero bias, for -Bj = .03 eY, using eqns. 

(4.46) and (4.48), and the values of and D for 
N D = 2.2 x 10- /cm 3 (given in Table 4.2), we see that in 
this case the condition that Q gc » Q f does not hold good. 

Since for the occurrence of T Q effect, we require 

^sc ^f the observation mentioned above provides a 

reason for non-observance of T q effect in cases where 

Eg, 0.03 eV for bhe values of other parameters as consi- 
dered here. 

To summarize the fir dings of this chapter, -e briefly 
mention the main results that have been obtained. 

(1) It has been shown that I Q effect can also be obtained 
- when the surface states are due to the tunneling of 
metal electrons into the forbidden bandgap of the 
semiconductor provided suitable values of and 
doping densities are chosen. 



* have succeeded la linking the characteristic 
energy S Q m the exponential surface state model with 
E t in the Belne-Pellegrini model, through a linear 
relationship between E Q and 1 /eJ. 



CHAPTER 5 


DETERMINATION OF DOPING- DENSITY FROM THE FORWARD 
' ' CURRENT-VOLTAGE CHARACTERISTICS 

In the previous chapter, the interest has been to 
understand the mechanisms of current transport that consti- 
tute T effect in Schottky barrier diodes. In this chapter 
we present a novel method for determining the impurity 
concentration in a semiconductor wafer. The method requires 
the measurement of forward I-V characteristics of a 
Schottky barrier diode fabricated on the semiconductor 
wafer mentioned above. One more requirement of this method 
is that the Schottky barrier should exhibit T Q effect. 

Note that this requirement automatically puts a limitation 
on the application of this method. 

Before presenting the methiod mentioned above, let us 
briefly go through some of the methods that are presently 
utilized to determine the impurity concentration of semi- 
conductor wafers. The existing methods can be broadly 
classified into two categories. The first category consists 
of those methods (e.g., differential capacitance versus 
voltage method) which can determine the impurity concentration 
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“ the as a functi “ Of depth [1.2]. Ihe second 

=ate g0 ^ consists of those Methods which can deterge the 
impurity concentration only when the doping density is 
independent of depth t2] . It is important to note, however, 
that the methods of category one are more versatile and can 
definitely do the 3 ob done hy the methods of category two 
fhe method evolved by us belongs to the second category 
and hence is useful only when the impurity concentration 
m the semiconductor wafer is independent of depth. 

1’tas we see that the method suggested in this chapter 
suffers from the following limitations: 

(1) It is applicable only when the Schotthy barrier diode 
-de on the wafers whose impurity concentration is to be 
measured exhibits a effect, and 

(2) It is applicable only when th e in- 
dent * doping density ie indepen- 

dent of -che depth. 

Hote that the method being suggested here utilises «. I-y 
characteristics. Xhe advantage of this method primarily 
ensues because of the ease of measurement of the forward 

I~v characteristics. 

Perhaps the most ■pooulan n-p +-u^> - j_ 

popular of the existing methods of 

category two is the measurement 

e measurement of capacitance -voltage (C-V) 

characteristic of Schottky barriers 

^ arriers on wafers whose 

impurity concentration is to be determine* „ 

Renamed. We consider this 

method m the next section. 
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5 . 1 E valuation of Doping Density from the C-V Characteristics 

3?he theoretical derivation leading to the use of C-Y 
measurements as a profiling technique is based upon the 
validity of the depletion approximation. This assumes that 
the depletion region consists entirely of charge due to the 
ionised donors, (i.e.) the mobile carrier concentration is 
negligibly small compared to the fixed impurity concentration. 

For a Schott ky barrier diode under depletion approxima- 
tion, the total charge Q gc on the semiconductor at any bias 
V is given by, 


Q sc = V b -V-4 f - TW/qJ! (5.1) 


where Up = the doping density 
and <i> B = barrier at a bias Y. 

Differentiating eqn. (5.1) with respect to the applied bias, 
we obtain the capacitance per unity area C' as, 


dQ 


kf> 


'V. 


- - m - + f > > 


(5.2) 


From eqn. (5.2) we observe that we require the voltage depen- 


dence of the barrier height to evaluate an expression for 
the capacitance as a function of the applied bias. Depending 
on the voltage dependence of the barrier height the diode 
can be classified either as an ideal diode with 3<J> B /3V = 0 
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or anon-ideal (practical) diode, whose non-ideality is 
defined by an ideality factor n, given by, 

n = [1 - 3< , b/)t] -1 ( 5 . 5 ) 

let ns now look into each of these cases individually. 

(l) Ideal diode: 

?or an ideal diode with _ 0 , fr0ffl eqn _ ^ 

ire get, 


O' = [ 


qe s II D V B 


( *B - Y ~*f * f 2 )]'* 


i . e. 


C ' 2 


qeN 


<*B - V ~ *f - s> 


-ll T-. 

S I) 


( 5 . 4 ) 

(5.4a) 


From eqn. (5.4a) we see that the doping density u can be 

Obtained from the slope 'm' of the l/C-’-V plot, usi ng the 
relation. 


D q* e 


s 


. 1 
1 0mm 

m 


(5.5) 


( 2 ) Non-ideal diode: 

" has been seen in Chapter I, the class of non-idea] 
can be classified into two categories depending on 
******* deB “dence of the ideality factor, fhese 

are, 

(i) modes showing Iq effect, wherein the ideality factor 
n Is given by the relation, 

n " (1 + V« (5.6) 
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where 2? Q is independent of temperature when measured under 
constant current condition. Note, at a given temperature 
is voltage dependent. 

(ii) Diodes showing *n* independent of temperature and bias 
'these diodes have a thin insulating layer between the metal 
and semiconductor and are referred to as MIS diodes.- 

Let us look into the C-V characteristic and hence the 
evaluation of N^ in each of these cases, 


(a) Diodes showing T q effect: 

For a diode showing r l‘ Q effect, using eqns. (5.6) and 
(5.3) in eqn. (5.2), we obtain, 

c ’ = c «f s E (,f B - Y - *± - r>]~ 7 < 1+ V T > c.7) 

From eqn. (4.13) in Chapter IV, we have, 

q 3 (J> td 

~ 2k(r|-l) ’ VlnQ (5.8a) 

SC 

q(4>^ - v) 

where p = - q4> f /kt (5.8b) 


Using eqn. (5.8a) and (5.8b) along with the relation. 


34>B 


E. 


FTnQ 


'sc 


in eqn. (5.7) and squaring we obtain, 


C 


I 2 


r 2 ( E ° , E o \i 


(5,8c) 


(5.9a) 



i . e. 


1 

O' 2 


1 


hr < x + it + 




(5.9b) 


/ ' kT 1 

where x = ($ B - V - $ f - |=0 

For reverse bias when E^/4xq 2 « (x + E 0 /q) we have [3]» 


and 


5* e s ^ N D * { n£ B " * f “ q 1 ) “ v ] > (5.10) 

3 4>b 


where * B = .4> BQ + — . Y = ^ + ^r 


Y 


flq = (1 + T o /T) 


T 0 = mean value of in the voltage range of 


interest. 


From eqn. (5.10) we observe that the doping density can 

_1 

C 

relation. 


be obtained from the slope , m t of the - V plot using the 


w = ~Z — . i • i~ 

E e g • q m n^ 


(5.11) 


(ii) MIS contacts: 

For these contacts the ideality factor 'n* is both 
temperature and bias independent [4]. fhus using eqn. (5.1) 
in eqn. (5.2) we obtain, 


1 

C 1 2 


3- e s^D 


* n [n(<{> B - 4 > f - ~) - V] 


(5.12) 


From eqn. (5.12) we observe that the doping density can be 
obtained from the slope »m* of the - Y plot using the 


C' 


relation, 
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N. 


D 


2 

1£ e 


1 

m 


n 


(5.15) 


Prom eqn. (5.11) and (5.15) we observe that to evaluate 

1 

the doping density Np from the ^ 5 - - V plot of a practical 
Schottky barrier diode requires a knowledge of the I~V 
characteristics also, Thus we see that it would be of 
interest to find a method to evaluate from the forward 
I-V characteristic and hence reauce the number of measure- 
ments to be carried out. 


Some of the other methods used to evaluate the doping 
density in a semiconductor wafer are (i) four point probe 
method and (ii) Van der Pauw method. Details of these 
methods are given in Appendix p* 

5.2 Determination of Doping Density from the Forward I-V 
Characteristics ~ 

Although forward I-V characteristics are convenient to 
measure and are extensively used to characterize a Schottky 
barrier diode, to the best of our knowledge, there has been 
no method available so far to determine the impurity concen- 
tration from the I-V measurements. The method to be des- 
cribed in this section is based on the dependence of T 
(calculated from the slope of the forward I-V characteristic^) 
on the current at which it is measured. For this we require 
a general relationship between T q and I in the case of 
diodes showing T effect. 
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Referring to eqn. (5.8a), one can write the expression 


for T o at a bias V as, 
2L = 


q 

2k( rj-l) * 31nQ, 


(5.13) 


'sc 


where 


q($ T,-V) 

“Td! (Nq/Rjj) 


q^ B -v) 


kT 


*- h = -[ln(l/SA*T 2 )+b] (5.13a) 


« doping density in the semiconductor 


❖ 


B = barrier height 


and 


*C ~ effective density of states in the conduction 


band of the semiconductor. 

Using eqns. (5.13), (5,13a) and (5.8c), one obtains. 


1 

T. 


2k 

«2L 


[ln(l/SA*T 2 ) + (b+l)] 


(5.14) 


O ' A ~o 

Prom eqn. (5.14) one observes that a plot of X/T versus 


ln(l/A*3? 2 S) will be a straight line of slope t - - and 


2k 


intercept with 1/T = 0 of, 


q»; 


Ini’ = -(b+l) + ln(A*T 2 S) 


i*e. , -b « (1 + In I f /A I 2 S) 


(5.15) 


(5.16) 


Using the relation b » In^/H^) » -la (U Q /N C ) and eqn. 
(5.16), the expression for is obtained as, 


% = N c exp (-b) = N (l»/A*T a S)-e 


(5.17) 


where e 


2.7183. 
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^ ^ Result s and Bisoussinn 

In Fig. ( 5 . 1 ) we present i- versus lnl plots at Z96°X 
for three Au-nSi Schottky barrie? diodes (BAu 110. EAu 120, 
EAU 130), which showed a ^ effeot (see seo _ 3>4> 

Ihe doping density obtained from the slope of the i- versus 

1111 Pl0t ’ USi “ e e «“- f 5 - 11 ) *« ^ese three diodes is given 
m fable 5.1, along with the doping density obtained from . 

the 0-7 method. The table also contains the values of doping 

density obtained for two more sets of Schottky barrier 

diodes. Bote that for the set of diodes made on bulk silicon 

" We “ 0lUded the <* ^ obtained from four point 

probe and Van der Pauw measurements also. 

Erom Table 5.1 it is observed that the value of doping 
density evaluated from the 1/T Q versus lnl plot matches well 
with the values obtained using the other methods. 

fo summarize, we have proposed a new method to deter- 
ge the doping density from the forward I-y characteristics 
ox a Schottky barrier diodes. With the proposal of a method 
evaluate the doping density from the forward current- 
voltage (I-T) characteristics, one observes that this method 

edong with the methods described in sec. 3.3 of Chapter III 

to evaluate E . a* n 

0 ^ B so» Can used to completely 

characterise a diode showing 3 > p-p-po,** ^ 

S i o efieot, from the forward I-V 




Plot of y versus in I for three small area diodes 
EAullO , E Au 120/ EAu130. 
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Table 5.1. Values of doping densities calculated using differ- 
ent methods for three sets of diodes showing f 
effect. 

loping Density (cm y ) 


probe } Pauw 


Wafer used 

Device 

1/T vs Ini 
°plot 

l/c 2 -V plot 

n + -n Sili- 
con Spj_ 
wafer 

EAu 110 

EAu 120 

EAU 130 

5.52 x 10 15 
5.35 x 10 15 
5.45 x 10 15 

5 3 io 15 
5.32X 10 15 
5.52x 10 15 

* - - 

+ 

n -n 
Silicon 

Epi wafer 

EAu 210 

EAU 220 

EAu 230 

EAu 240 

EAu 250 

5.0 x 10 15 

5.1 HO 15 
5.1 x 10 15 
5.08 x 10 15 
5.16 x 10 15 

5.2 x 10 15 

5.l6x 10 15 

5.05x 10 15 

5.17x 10 15 

15 

5.4 x 10 ? 

Bulk 

n-Si 

wafer 

DAu 111 

J>Au 121 

DAu 131 

DAu 141 

2.35 x 10 15 
2.25 x 10 15 
2.40 x xO 15 
2.30 x 10 15 

2.22x 10 15 
2.25X.10 15 
2.3 x 10 15 
2.1 x 10 15 

■ ■■■■ ■ ■ t 


2.4xl0 15 2.3xl0 15 
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data alone. This means a reduction in the number of 
measurements to be carried out on the diode. Note, 

-^so re P resen t respectively, the characteristic energy, 
the neutral level at the surface and the density of surface 
states at the neutral level in an equivalent Levine's 
model [5] for a diode showing T anomaly. 

Before concluding, we would like to reiterate that 
this method has the following limitations: 

(i) It can only be used if the Schottky barrier diode 
shows T q anomaly, and 

(ii) It can only be used if the semiconductor wafer has a 
uniform doping profile. 



CHAPTER 6 


SUMMARY AMD CONCLUSIONS 

It has "been, observed by many workers that the ideality- 
factor (n) of rectifying metal-(low doped) silicon ( or CaAs 
etc.) contacts can either be independent of or dependent on 
temperature (!') depending upon the processing and the nature 
of starting materials. The work in this thesis has mainly 
been concerned with those diodes which show a temperature 
dependent ideality factor. The main results are summarised 
below. 

(i) We observed that a method has been suggested 
in literature to distinguish between the thermionic emission 
(TE), Thermionic field emission (TEE), and field emission 
(PE) as the mechanism of current transport in Schottky 
barrier diodes. The method consists of plotting nT versus T 
(obtained from the forward current-voltage (I-V) characteris- 
tics measured at different temperatures) where a linear plot 
corresponds to TE, a temperature independent plot corresponds 
to PE and the intermediate case corresponds to TPE. The 
work done here has pointed out that the above mentioned 
assignment is valid only in the case of homogeneous diodes 
(where there is no fluctuation in geometrical and material 
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parameters transverse to the metal-semiconductor interface) 
and is not valid for nonhomo gene ous diodes. For example, 
it has been shown, both theoretically and experimentally 
(see Chapters I, II and III) that if T q (where I is a 
temperature independent constant in the relation nI=I+I o ) 
fluctuates the nT versus 1 plot would indicate the current 
transport to be due to field emission in the conventional 
sense when it is actually due to thermionic emission. I'hus 
the present work emphasizes that caution must be exercised 
when inferring the mechanism of current transport from the 
nature of the nT versus T plots. 

(ii) in attempt was made to correlate f Q with the 
material and geometrical parameters at the metal -semiconductor 
interface. In this regard it was noted that Levine (in 1971) 
had already proposed a model where an exponential energy 
distribution of interface states, characterized by a 
characteristic energy E Q? was shown to be the reason for 
the appearance of T q effect in Schottky barrier diodes. 

Later in 1977, Crowell has pointed out that the choice of 
exponential energy distribution of interface states is not 
the only one which would give T Q effect, rather several 
other distributions of interface states could do it. For 
example, he has shown that it is possible to obtain T effect 
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in M-P( thin) -H diodes, provided suitable doping profile 
is chosen in the P~layer. We further noted that Heine 
(in 1965) has proposed a model of Schottky barrier where a 
potential profile similar to that of M-P( thin) -N diode 
appears due to the charge incorporated by metal electrons 
tunneling into the forbidden bandgap of the semiconductor. 
Prompted by these observations and the detailed calculations 
based on Heine's model as carried out by Pellegrini (in 1974), 
we set out to correlate a parameter E,!,, the limiting value 
of the transverse energy of the electrons tunneling from 
the metal into the forbidden bandgap of the semiconductor, 
with the T q effect. The calculations have been done in 
two steps. She first step was to assume the correlation 
of T q with E Q as given by Levine, and the second step was 
to correlate the intermediate phenomenological parameter 
E q , with another phenomenological parameter El. It has 
been observed, as shown in Chapter IV, that Ey and E Q can 
be correlated provided and the doping density of the 
substrate lie within a specified range. Thus, we note, 
that fluctuations which modify the n!i versus I plots in 
the case of thermionic emission, so that apparently field 
emission seems to be the mechanism, may be due to fluc- 
tuations in E^,. Phis provides at least one physical 
mechanism to explain the fluctuations in T . 
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The absence of the identification of a physical 
■mechanism to explain the exponential energy distribution of 
surface states, and the identification of T Q effect in 
cases other than exponential surface states leads one to 
the question, whether Levine's model is only a mathematical 
representation of the actual conditions existing at the 
interface. A definite answer to this question would require 
an extensive study of the metal semiconductor interface 
and the formation of surface states. Accepting Heine- 
Pellegrini model to be the source of T q effect, we would 
require the completion of the following studies to identify 
the cause of fluctuations in E£ and hence T 0 : 

(a) an experimental verification of the existence of T Q 
effect in Heine -Pellegrini model, 

(b) a theoretical evaluation or an independent experimental 
determination of E^, and 

( c) the modifications required in the He in e-Pellegrini 

model in the presence of a thin oxide layer between the 
metal and the semiconductor. 

(iii) A method has been developed to determine the 
doping density of a semiconductor wafer using the dependence 
of T on the current at which it is measured. This method 
requires the measurement of only the forward I-V charac- 
teristics of a Schottky barrier diode, hence can be used 
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as a convenient method, if one does not want to use the 
conventional method of determination of doping profile 
from C-V measurements. The limitations of the method 
based on forward I—V characteristics are given in C ha pter 5. 

This bring us to the end of the present work, which 
was mainly concerned with the understanding of the mechanism 
of current transport in nonhomo gene o us , rectifying metal- 
semiconductor junctions showing temperature dependent 
ideality factor. 
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APPENDIX A 

EVALUATION OP 4> bo and n FOR AN MIS DIODE 

In this appendix, expressions are derived for the zero 
bias barrier height <J> B q and the ideality factor n for a 
metal-semiconductor contact with an interfacial oxide layer 
of thickness 5 and a constant energy distribution of surface 
states. Pig. (A.l) shows the energy band diagram for such a 
contact under (a) forward bias, (b) zero bias, and ( c) reverse 
bias conditions. 

The surface states are considered to be of two kinds, 
one controlled by the metal Permi-level , where the surface 
state charge is given by, 

B 

0 = -q 2 D (-£_$__$) ( a-1) 

ssa H sa v q B 0 J ' 

where D ga = density of surface states controlled by 

metal -Fermi level in cm eV . 
q$Q = neutral level measured from the top of the 
valence band in eV. 

• q.4>g = barrier height in eV 

and Eg = bandgap in eV. 

and the other controlled by the semiconductor Permi-level 

where the surface state charge is given by, 

= -'jXb <i* - *B - *0 + V > 


( A-2) 





176 


where V is the applied bias with the proper sign. The total 
surface state charge is then given by, 

%b = <W + Q S sb “ <5® * % - V 

E 

+ q2i) sb (A ” 3) 

The space charge in the depletion layer of the semiconductor 
is given by, 

Q sc = V2«ie s N D T* B - V - TWOT (A-4) 

The. total surface charge density, on the semiconductor surface 
is given by the sum of eqn. (A-3) and (a- 4) and for charge 
neutrality in the system we can write, 

= -<«es + »sc> ‘ 

where is the charge on the metal. 

Considering a potential drop A across the interfacial 
layer, and the application of Gauss's law to the surface 
charge on the metal and semiconductor, the charge on the metal 
is obtained as, 

-e.A 

«» = ~r~ < A - 6 > 

where is the permittivity of the interfacial layer. From 
an inspection of the energy band diagram of Fig. (A-l) we 
have, v 
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4 (t> m - X ox ) + 94 = l(x s - X ox ) - 'lv g 

i.e. 4 = <t m -X s - (A-7) 

Using eqns. (A-3), (A-4), (a- 6) and (a- 7) in eqn. (a- 5) , we 
obtain. 


(4> 


/ 24E s N d 62 


131 


x s - V 


<*B “ V “ *f - f ] 


(q^ " " V + ^ 2j£) s b “ * B “ * 0 + 

( A-8) 

The barrier height <i> at zero bias can be obtained from 
eqn. (A-8) by putting V=0 and solving for <p B . The solution 
for eqn. (A~8) with V=0 is obtained as. 




E. 


BO 


C°2 <♦ m "V + (1 - C 2 ) (-* -♦„)] 


°I °1 


+ *** - ^ /a t%(*-- ,l ») +( l- 0 2)(S K -*o) ^ 

°1 , kT °o C i i 

- ^ (+f + + -pk' 


( A-9) 


where C-, 


2 W 


5 o 2 = 


e i 


< e i + ^ 63 > ss > 


and I) = B + B . . 
ss sa sb 


Our interest now is to obtain the ideality factor n, 
which is given by the relation. 



l 


n 


■ 


(A-10) 

Differentiating eqn. (A-8) with respect to V, we obtain, 


5 ^B _ 6 r e s , 3 *B 


34> 


[ ¥‘ h T - X) + 4 “sa *T 


D_ 


£ + i 2;u S b <ir - ^ a- 11 ) 


where ¥ = depletion width = 


2e 


s 


. c i ir D (4> B 


Y - 




f q. 


Adding 1 to both sides of eqn. (A~ll) we have. 


9<J> -3$,, e q 2 B 6 

- rfi = It - i) + ***„ + «\b> + -rf- + 1 


T> 

*.* (l - g - y ~) = — , eqn. (A-12) reduces to 


1 

n 


1 + q 2 IL 6/e. 
_ sa j x 


(%r + q2;D -~ + I- + i 


sa 


sb e i 


(A-12) 


i . e. n = 1 + 


(e 3 /V t S 2 D ab > 6/e ± 


(1 + l^sa 6 /e 1 ) 

Bor e g /¥ <<: q 2 B s ^, eqn. (A-13) reduces to 


(A-13) 


n = 1 + 


q z D 


sb ‘ 


(1 + l! “ aa 6 / £ i ) 


1 + ^sb 6 / e : 


( A-14) 


where 


sb 


(1 + 9 2 B sa 6/0 


Thus we see that eqn. (a- 9) and (A-14) give us the relation- 
ship fortfcgQ and n in the presence of two kinds of surface 
states for an MIS diode. 
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APPENDIX B 

AN ALGORITHM TO FIND THE SMALLEST ROOT OP A POLYNOMIAL 


In this Appendix we present a method to evaluate the 
smallest root of a polynomial correct to two decimal places. 
The root finding algorithm is based upon the following 
theorem that was initially derived to yield an approximate 
root value in terms of all the coefficients of an arbitrary 
nth order polynomial (L. Tessler and L. Eisenberg, Proc. IEEE, 
Vol. 60, p 737, 1972). 

Theorem: Given an nth order polynomial 

PjjCx) = 1 + AjX + a 2 x 2 + ... + A^ ( B-l) 


the smallest magnitude root of (G-l) has an approximate 
value. 


R 


= ~1/(A 1 + JJj) 


( B— 2) 


where 


h = 


-a 2 a 


0 

1 


0 

0 


“% %-2 


0 

0 


0 

1 

An 


¥ 


-1 


( B-3) 
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A^ 1 0 0 .... 0 

"^2 A 0 0 

• 

! o 

. i 

%-l ' A N-2 A 1 


{ B-4) 


Since the elements in eqns. (B-3) and (B-4) appear in a 
systematic manner, it is possible to construct an algorithm 
with very few arithmatic operations. Specifically if we 
define , 


A ± 1 0 

a 2 A^ 1 0 





(B-5) 


and Wq = 1 

then Wjj can be calculated from the recursion relation, 

% = (-1) 1 (B-6) 

Using eqns. (B-5), (B-3) and (B-4) in (B-2) one gets, 

E = - l/< VV = - (B-7) 
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In practice, we see that most of the polynomials are of 
the form 

Q(x) .= C N x N + C N x 5 " 1 + .... + C 0 , C 0 ^l (B-8) 
For polynomials of this kind the determinant counterpart of 
% ls 

0 1 0 0 0 0 

^2 ^1 ^0 • # • * * Q 

♦ 

* 

C N -1 * °0 

C* N w 

°o % 

which reduces eqn. (B-6) to, 







( B-9) 


It has been shown by 1. lessler and 1. Eisenberg (Proc. IEEE, 
Vol. 60, p 737, 197 2) that eqn. (B-9) gives the value of the 
smallest root correct to two decimals. Since we are not 


interested in very accurate values, we have chosen this 
algorithm to evaluate the root close to the origin in the 
cubic equation encountered in Chapter II. 
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APPENDIX C 

SOLUTION OP TWO DIMENSIONAL POISSON* S EQUATION EOR AN 
AUXILIARY SCHOTTKY BiRRIER GUARD RING STRUCTURE 

In this Appendix we give the details about the choice 
of the boundary and the boundary conditions for the solution 
of two dimensional Poisson’s equation in a Schottky barrier 
diode with an auxiliary schottky barrier guard ring structure. 
The presentation is divided into two parts, viz. (a) the 
formulation of the problem and (b) the solution of the 
Poisson's equation. 

(a) Eormulation of the problem: 

The Poisson's equation in two dimensions, considering 
the mobile carriers (majority carriers only) and the charge 
due to the ionised donors (N^) can be written as, 

+0 - -^[l-esp(-»/M>] (0-1) 

3X 2 SY 2 e S 

where $ is the potential at any point described by the co- 
ordinates X and Y. The next step now is to identify the 
boundary and the corresponding boundary conditions. Prom 
a perusal of Pig. (C.l), where we present a cross-sectional 
evaluation of a Schottky barrier diode, we see that the 
lines pp’ , pr, ps get automatically chosen as the boundaries, 
whose pp* forms the line of symmetry, the line pr defines the 




Cross-sectional elevation of a Schottky barrier 
diode with an auxiliary Schottky barrier 
guard .ring. 


Fig.C.2 Schematic drawing of the mesh used in 
the numerical computation. 


i» 
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line of metal-semiconductor junction and ps defines the 
hack contact. The fourth boundary is chosen as follows. 

Let us consider the line qq' perpendicular to the junc- 
tion and passing through the midpoint between the main and 
the auxiliary Schott ky barrier diodes. It has been shown by 
McKenna and Itfasserstrom (B.S.T.J. , Vol. 49, P 853, 1970) 
that a metal strip of half width greater than thrice the 
depletion width (W) under the metal can be considered to 
be semiinfinite for potential calculations. Based on this, 
qq' has been considered to be the fourth boundary (as a line 
of symmetry) assuming the radius of the main diode and the 
half width of the guard ring to be greater than 3W. Thus, 
we note that the region for the solution of the Poisson's 
equation is pqq'p 1 in Pig. (C.l), and the boundary conditions 
can be written as follows (see Pig. C.l) 

(i) Along pq (X=0) 

$(0,Y) = (<f> B - V - 4> f ) =$ 0 0^ Y^ r (at the metal) (C.2a) 

and using the criterion of continuity of potential and the 
normal component of the displacement vector 

»»(o + .Y) = (Q-.Y) 

T) ITj [ 3 X 

r< 14 (r+d/2) (C.2b) 

®(0 + ,Y) = 


* (0“,Y) 
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(ii) Along qq* (from symmetry condition) 

= o (c-3) 

(iii) Along q' p* (back contact) 

*(t,T) = 0 ( C-4) 

(iv) Along pp' (symmetry condition) 

3£inol . 0 (0-5) 

for the above mentioned mixed boundary condition, the 

two dimensional nonlinear Poisson's equation given in ( C-l) 
does not have an analytical solution and hence Numerical 

methods were resorted to for the solution of eqn. (C-l). 

• ; ■ .. , ■ y; ’ , , Vi v 

( b) Solution of the Poisson's equation: 

Por easier handling, let us first normalize the dis- 
tances with respect to the depletion width ¥, the potential 
with respect to the band bending (^ - V and the 

charge with respect to qN-p, 

where ¥ = 2Xp VJ^T 

% = # 0 /OdJ/q). (C-6) 

and = I)ebye length. 

Carrying out the normalization mentioned above eqn. (C-l) 
reduces to , 


j 
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a 2 <fr + ill = 2[1 - exp(-l$ a l$)] = p 

3x 2 ay 2 

The eqn. (C-7) is now numerically solved by replacing 
the differential equation by a five point difference equa- 
tion and solving the resulting equation by the Newton-Gauss 
implicit scheme with a step size of h = 0.05 and the 
relaxation factor w = 1.73. 1'he differential equation then 

reduces to. 


^l+l,3 + *i-l,j +4 'i.3 +1+<l ’ i '3- 1 ' 


= - f- [1 - exp( ~ l*el * *ij^ 

leading to an iterative relation, 

^(n+i) _ |i exp (- |4> s I ♦ 


(C-8) 


- £ [♦■<*> • + * + *H+i + ♦ - 2h2](c ~ 9) 

“ 4 L *i+l,j y i-l» 3 1 »3 + - L 


Let us now look into the boundary conditions and hence the 
iterative relation at the boundaries. Fig. (D.2) represents 
the given structure divided into mesh of dimension MXN with 

the mesh size equal to h. 

(i) Since the two dimensional nature of the potential 
profile is restricted to within three times the depletion 
width from the edge of the metal (B.S.T. J. , Vol. 49, P 853, 
1970), beyond the line P* P we have, 
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\ = 0 (C- 10 ) 

Based on this observation we shift the boundary at pp‘ to 
PP* to reduce the number of mesh points* The iterative 
relation along the line PP' (i=0) can then be written as, 

<K0,j) = $(2,3) 1 ^ ^ H (c-11) 

(2) Along the line qq* (i=M) we can write 

<fr(MHL, 3 ) = <J> (M-l, 3) 1^ 3$ N ( 0 - 12 ) 

(3) Along the line Pt, where the front metal is, we have, 

<{>(i,l) = 1.0 14 ±4 Ml (0-13) 

The boundary condition along the line tq (see Fig. C.l) 
is in a little unwieldy form to be applied for the numerical 
calculations, for it would need a very small mesh size near 
the boundary and hence would lead to an unnecessarily large 
number of mesh points. To get over this one replaces the 
boundary condition (C-2b) by its equivalent integral form, 

j e dl = / P dxdy (0-14) 

G 3 n # G 

where G is the path chosen by connecting the midpoints of 

the surrounding mesh (see Fig. C.2). This reduces numerically 

to, 

<J> (i,2) + r)4>( i , 2) + ($(i+l,l) + 4> (i-1,1)) - 2(l+r)) 4>(i»l) 

= h 2 [1 - exp(-}<j> a | <j> (i,l))] (0-15) 
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Using this in the iterative relation we obtain 
^ - h2 


i rA^ n+1) + + A ( n ) i~ 

4 L< ** i-1,1 $ i+1,1 $ i,2 J 


( 0 - 0 . 6 ) 


(4) The fourth boundary at the bach contact is shifted to 
a distance ofstimes the depletion width from the metal- 
semiconductor junction. This has been done to reduce the 
number of mesh points. The boundary condition used is 


,,N. “ $ i , F+l 


M 


(C-17) 


and j=N represents the boundary. 

With the boundary conditions outlined above, the Poisson* s 
equation is solved using the scheme, 


4, ( n+1 > m £ ( n \ - M 

w i,j. v i,3 4P 


where F 


ih.j + ♦ ' 7, j+i 


(C-18) 

(n-1) (n+1) M (n) 

+ <? 4 _/l 4 + 4 > -i 4_1 “ w 4 4 


- 2h 2 [l - exp C- $ If] 3)1^0 . 3 


(C-19) 


= -~Tzr~ = -I** 2 N> s i^ exp 1,3 !♦.(> I < G - 20 > 

9< J> i, j 

w = Relaxation factor = 1.75. 


and the convergence criterion, 
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or 

The 

and 


6 i, 3 


— W ( £? - -* l t i i < 10 ' 4 


6, , = /WV ) 2 /(MM+M) < 10~ 5 


mesh size chosen for our work was M = 101, Ml = 61 
U=6l. 


(C-21) 
( 0 - 22 ) 
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APPENDIX D 

EVALUATION OP THE POTENTIAL PROFILE IN THE 
HEINE- PELLEGRINI MODEL 


In this Appendix, we present the method for evaluating 
the potential at the metal -semi conduct or interface and hence 
the complete potential profile in the Heine -Pellegrini 
model. As seen in Chapter IV the Poisson’s equation in a 
Schottky barrier diode can be written as, 


[N v exp(-x/X v ) + N c exp (-x/X c ) - N^] (D-l) 

CLX. q 

In (D-l) the densities of metal electron N B tunneling close 
to the band B of the semiconductor and X the corresponding 
penetration depth are given by (here B=c for conduction band 
B^v for valence band) , 


N B = 


E osB E T ^xB ^TB ln ( 1+A B2 - 


2E m 


372" 


1+A- 


B1 


ft {arc sin ( . - - tt) - arc sin[ +1 

A B2 +x a B] _ 


( 1+1 /A B2 ) A B1 -1 


X B “ 


■3} 


2 ^ m xsB tA xsB S osB ( 1+ ^A B2 )] T ln(l/A B2 ) /(l+A-^) ] 


Aqi m ~ *s “ V + ^T^^^^xc^osc^ 


m. 


A, 


xsV 


B, 


T 


C2 


' m TT — m 
xs V xs c 


) /(u E ) 
;/x ^XC oscr 


(D-2) 


(P-3) 


(D-4) 


in the conduction band and, 



A yl - 0.0 

m E Em • 

^2 - ( I^§5g - T + #> I <1 t'xvl W < B -5) 

In equations (B-2) to (B-5) the symbols have the fol- 
lowing meaning: h = 27th is the Plank constant, m^-^ is the 
electron effective mass in the band B of the material H 
(H s m and H = s in the metal and the semiconductor respec- 
tively) in the direction D (in obtaining (B-2) and (D-3) 
it is assumed that the constant energy surfaces are ellip- 
soids with the revolution diameter coinciding the x axis so 
that B = Z represents x direction and D = 1 indicates a 
direction perpendicular to the x direction), 

Ejyj is the metal Permi -level evaluated from the 'bottom of the 
metal conduction band, is the metal workfunction, X g is 
the semiconductor electron affinity, E is the amplitude of 
the semiconductor forbidden energy gap, B osc = B c ** 

+q(4» m - X g ) and E ogT = E v + E c - E^, are the electron 
energies at edges of the forbidden gap, v Q is the potential 
at the metal-semiconductor interface and E^ is the limit 
value of the electron energy perpendicular to the x axis 
beyond which the quantum penetration into the forbidden 
energy bandgap of the semiconductor is negligible 
(B. Pellegrini, Solid State Electron., Yol. 17, P 217, 1974). 

A closed form solution for eqn. ( D— 1) has been obtained 
by assuming that there exists a length w of a section of the 



192 


depletion region characterised by the Inequalities, 
v - V - v(w) >> kl/q 

(D— 6) 

and W » w »X Q , X^. 

where v c = 4 m ~x 3 ) - is the contact potential difference 
V = applied bias 

v(w) = the potential at a distance w from the metal- 
seini conduct or interface 
and ¥ = depletion width. 

Under this assumption the field "S(x) , and the potential 
v(x) at any point within the semiconductor, at a distance x 
from the metal -semiconductor interface can be written as, 

t-(x) p;M [I c X c expC-X/X ^} + N V X V exp(~x/x v ) - K^V-x)] 

(D-7) 

q¥f,x ■; , 2 

V(x) = v 0 — (¥ - |) + §- [1 - exp(-X/^ C )]N C X Q 

>3 s 

+ [1 _ exp(-x/X v )] K v x ^ (li-8) 


where 




the contribution of the surface state charges 
to the built in potential 


n [h c xJ + 


(h-9) 
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v e = the contribution of the electron tunneling from 
the conduction band into the forbidden bandgap 
kT 

Our interest now is to obtain the value pf v 0 and hence 
the potential profile, and is done as given below. 

Invoking the condition of continuity of electric flux 
at the interface, we have, 

e m c(0") = e s t(O^) (D-10) 

where & (0~) is the field in the metal for x=0~ and f(0 + > is 

*F 

the field in the semiconductor for x=0 * If one neglects the 
positive charge induced by the band bending, the field75(0~) 
in the metal can be written as, 

'S(o-) i + Jjl| ( B— 11} 

IF 

Using (l)-ll) in ( D-IO) , we obtain the relation, 

= q [H c X c + h v x t - X D V ] (D— 12) 

It has been shown by Pellegrini ( J. Phys . I): Appl. Phys., 

Vol. 19, P 55, 1976) that in most of the metal-semiconductor 
systems, N » 0, this redupes eqn. (D-12) to, 

V 

- qUpW) (h-13) 


e v 
m o 


TI 
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Since ¥ is a function of v Q , eqn. (0-13) is a transcendental 
equation and has been solved numerically by using iterative 
methods . 

After having evaluated v Q , our next step has been to 

evaluate the potential profile and hence the barrier height 

b in the cases considered by us we have£(0 + ) > 0, thus 
B° 

leading to a maxima in the v(x) profile. The maximum value 
of v(x) = v m occurs at a distance x m from the metal-semi- 
conductor interface. The barrier height <*> B is then given 


by. 


4>B = <♦« " X s> + v 


m S' 


m 


(0-14) 


The expression for x m and hence v ffl can be obtained as follows. 

At x = x , we have (x) =0, hence from eqn. (0-7) for N c = 0 
m 

we have, 

"€(* m ) = £- IV, *«p - n u (w-^ m )] = o (D-i5) 


"S 


i.e. exp (-x m /X T ) = (1 - x^/V) 


( 0 — 1 6 ) 


Taking natural logarithm on both sides series expanding 
ln(l-x m /W) for x m /W « 1 we have, 


m (l»A») = - vT + 


v /M D 

x 'w 


v 


i-e- = (x',-VT ln »v vV 0 


(D-17) 



Using eqn. (U-17) in eqn. (.0-8) the potential at 2 = x m 
and hence the barrier height can be obtained. 

In our calculations we have considered E m , the metal 
Perm! Level evaluated fro:, the bottom of the metal conduction 
band, and N , the density of electrons in the metal conduc- 
tion band, to be temperature independent. The temperature 
dependence of 4> m , the metal workfunction, X g , the electron 
affinity of the semiconductor and E , the bandgap of the 

O 

semiconductor are given in Chapter II. 
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RESISTIVITY MEASUREMENT BY POUR POINT IROBE 
INI) VAN PER PAU¥ METHOD 


the Van der Pauw methods to determine the resistivity of a 


semiconductor -wafer 


Pour point probe method: 

A schematic of the dc circuit used with a four point 


probe is shown in Pig 


% **>.*'[' , ; 


Constant Current 
Source 


Semiconductor 


Schematic showing the dc circuit used with a 
four-point probe method. 
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por a probe resting on a semiconductor with the nearest edge 
at a distance more than 4s from any probe to the nearest 
edge, the resistivity is given by, 

(E-l) 


P= 2hb • X ' *(*/«) 


In eqn. (E-l) a represents the probe spacing and t 

, on A F is a correction factor which 

thickness of the samp » 

, the (t/s) ratio . Normally the probe spacing « 

depends on the W 8 ' thick- 

hnnt 16 cm, so that 2its = 1, an 
chosen to be about; .k > 

ness t > 33, ean. (*>-1) redU ° eS t0 

y (*■*> 

p - y ohm-cm 



L-CUi- 


A He mentioned here that care mu» . 


[ surface leahage and that g 

made 


StHfp 

PI 


(b) Van der Pauw method : 

1 di'VXji, 


th dT measure the resistivity of a semiconductor 
A method to measu 


A. ffit# uuuu , , 

. . nf its size and shape was proposed by 
wafer , irrespective of its size 

■ iorr (Philips Research Report, Vol. 13, 

Van der Pauw in 1958 (mini* 
p 1, 1958) and is as follows. 

p ls E.2 shows a sample of arbitrary shape used with the 
van der Pauw method. *» -pie has four small contacts a 
the circumference of the sample. The resistivi y P 
ss nmole is given by, 


of the 



Fig. B.2. 


A sample of ‘arbitrary shape used the Van der 

Pauw method. 





In eqn. (E-4a) Vqj> 
G and D for a cur r 


through the 


that they are syumewx^ 
tadjacent contacts, 

>s function also hecomes 




equal to 1. It should, be noted that the sample for this 
method should be of uniform thickness. 

Before concluding, we would like to mention that the 
methods mentioned in this Appendix are not applicable to 
n/n + or p/p + epitaxial wafers, since in these measurements 
the epitaxial layer gets shorted out by the heavily doped 
substrate. 



